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ABSTRACT 

The access of 1.2-40 MeV protons and 0. 4-1.0 MeV electrons from 
interplanetary space to the polar cap regions has been investigated with 
an experiment on board a low altitude, polar orbiting satellite (0G0-4). 

A total of 333 quiet time observations of the electron polar cap 
boundary give a mapping of the boundary between open and closed geomag- 
netic field lines which is an order of magnitude more comprehensive than 
previously available. 

Persistent features (north/south asymnetries) in the polar cap pro- 
ton flux, which are established as normal during solar proton events, 
are shown to be associated with different flux levels on open geomagnetic 
field lines than on closed field lines. The pole in which these persis- 
tent features are observed is strongly correlated to the sector structure 
of the interplanetary magnetic field and uncorrelated to the north/south 
component of this field. The features were observed in the north (south) 
pole during a negative (positive) sector 91% of the time, while the 
solar field had a southward component only 54% of the time. In addition, 
changes in the north/south component have no observable effect on the 
persistent features. 

Observations of events associated with co-rotating regions of en- 
hanced proton flux in interplanetary space are used to establish the 
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characteristics of the 1.2-40 MeV proton access windows: the access win- 
dow for low polar latitudes is near the earth, that for one high polar 
latitude region is ^250 behind the earth, while that for the other 
high polar latitude region is M750 behind the earth. All of the 
access windows are of approximately the same extent (M20 R ffi ). The fol- 
lowing phenomena contribute to persistent polar cap features: limited 
interplanetary regions of enhanced flux propagating past the earth, 
radial gradients in the interplanetary flux, and anisotropies in the 
interplanetary flux. 

These results are compared to the particle access predictions of 
the distant geomagnetic tail configurations proposed by Michel and Dessler, 
Dungey, and Frank. The data are consistent with neither the model of 
Michel and Dessler nor that of Dungey. The model of Frank can yield a 
consistent access window configuration provided the following constraints 
are satisfied: the merging rate for open field lines at one polar neu- 
tral point must be ^5 times that at the other polar neutral point, re- 
lated to the solar magnetic field configuration in a consistent fashion, 
the migration time for open field lines to move across the polar cap 
region must be the same in both poles, and the open field line merging 
rate at one of the polar neutral points must be at least as large as 
that required for almost all the open field lines to have merged in 
0{or\e hour). The possibility of satisfying these constraints is investi- 
gated in some detail. 


The role played by interplanetary anisotropies in the observation 
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of persistent polar cap features is discussed. Special emphasis is given 
to the problem of non-adiabatic particle entry through regions where the 
magnetic field is changing direction. The degree to which such particle 
entry can be assumed to be nearly adiabatic is related to the particle 
rigidity, the angle through which the field turns, and the rate at which 
the field changes direction; this relationship is established for the 
case of polar cap observations. 
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I. INTRODUCTION 


An investigation of the physical processes involved in the access 
of low energy particles from interplanetary space into the interior of 
the earth's magnetic field is of special significance to several geo- 
physical and astrophysical questions. Such an investigation bears 
directly on the question of the configuration of the distant geomagnetic 
field and the relation between this field and the polar cap region. A 
definition of this configuration is perhaps the only means available for 
the determination of the degree to which magnetic field merging plays a 
rble in the interaction between solar and terrestrial plasmas and mag- 
netic fields. Constraints placed on the extent of magnetic merging in 
this interaction would have far-reaching implications for other phenomena 
involving solar magnetic fields and for some phenomena involving galactic 
magnetic fields. 

Since the discovery, in 1964 [1], that the solar wind distorts the 
geomagnetic field to the extent of forming a "tail" in the anti -solar 
direction, there has been a good deal of speculation concerning the con- 
figuration of the magnetic field in this tail at large distances from 

I 

the earth. The models which have been presented for this configuration 
by Dungey (1961) [2], Dessler (1964) [3], and Frank (1971) [4] differ in 
their estimates of the length of the tail by one or two orders of magni- 
tude. As Dessler pointed out, charged particle observations in the polar 
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caps constitute one of the most efficient and appropriate means of deter- 
mining the structure of the distant geomagnetic field, since it is not 
practical to use satellite-born magnetometers to map this structure in 
detail. 

Since the fundamental theoretical difference between the three 
models of the distant geomagnetic field is the assumption made about the 
extent and/or mechanism by which the solar and terrestrial fields merge, 
a resolution of questions concerning the configuration of the distant 
field can have a direct bearing, through these models, on the question 
of whether magnetic field merging is a significant process in the inter- 
action between the solar wind and the geomagnetic field. Since it has 
been suggested that magnetic field merging may play an important role in, 
for instance, the generation of solar [5] and/or galactic [6] flares, a 
determination of the possibility or impossibility of magnetic merging in 
astrophysical plasmas, even over a limited range of plasma parameters, 
would be quite significant. The interaction between the solar wind 
plasma and the geomagnetic field is the only readily available system 
for which this can be determined. 

Although charged particle measurements cannot determine the extent 
of magnetic merging directly, they can lead to the establishment of con- 
straints to be placed on magnetospheric models. Several studies have I 
been conducted with this goal in mind. Observations of >50 keV and 
>20 keV electrons in the magnetotail have been interpreted by Van Allen 
[7,8] and Anderson and Lin [9] as evidence that these electrons gain 
access to the magnetotail along geomagnetic field lines which are 
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connected to the interplanetary field. These results have, however, been 
interpreted differently by Michel and Dessler [10], who have modified the 
closed field configuration model (e.g. no merging) to account for these 
electron observations. In any case, the mode of access for >1 MeV pro- 
tons may be quite different than that for >50 keV electrons, owing to the 
large differences in magnetic rigidity (MD.05 MV for the electrons vs. 

^1 00 MV for the protons) and velocity (0.43c (electrons) vs. 0.046c (pro- 
tons)). The observation by Evans and Stone [11] of a large north/south 
difference in the proton polar cap flux lasting for more than twenty 
hours showed that the question of proton access was still unresolved. 

Prior to the observations reported by Evans and Stone [11], some 
evidence of structure in observed polar cap proton fluxes was available 
[12-21], but the observations were too sparse to lead to comprehensive 
analysis. Since this preliminary report, several observations of per- 
sistent north/south asymmetries have been reported [22-26], and several 
of the most recent (Englemann, et at. [24], Van Allen, et at. [25], and 
Morfill and Quenby [26]) have been interpreted in the context of the 
Dungey open field configuration, relating the north/south polar cap flux 
differences to interplanetary flux anisotropies. Van Allen, et at. [25], 
for instance, report observations which seem to follow this relationship 
closely, but only for a period of about six hours during one solar flare 
event. This emphasizes a limitation which is common to all the previous' 
observations cited above: none of these studies deal with observations 
from more than two events. Results obtained from the analysis of only 
one or two flare events are subject to the severe limitation that the 
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solar and interplanetary parameters which determine the behavior of the 
particle flux near the earth can vary significantly from one solar flare 
event to the next. Thus, although it is now well established that large 
scale structure in the polar cap flux of low energy protons is not un- 
common, aside from the indications that interplanetary flux anisotropies 
play a role in the observations of north/south asymmetries, the relation- 
ship between polar cap fluxes and interplanetary particle fluxes is still 
unresolved. 

The high time, energy, and flux resolution of the data presented 
in this thesis make them particularly appropriate to an investigation of 
low energy charged particle access to the polar caps. As suggested by 
Vampola [27], low rigidity electron observations can be used to map the 
geomagnetic tail field lines onto the polar caps in order to determine 
the polar cap boundary between open and closed geomagnetic field lines; 
the 333 observations of this boundary reported here will result in a 
much more comprehensive mapping than the 25 such observations previously 
available [27]. The proton observations reported herein differ from the 
observations cited above in at least two important aspects: (1) the avail- 
ability of data throughout an eighteen month period has resulted in the 
compilation and analysis of fifty-four solar proton events rather than 
one or two, and (2) some of these events represent the first polar cap 
observations of fluxes associated with regions of enhanced flux in 
interplanetary space which are co-rotating past the earth. These obser- 
vations will be interpreted in terms of the locations of the access 
windows for 1.2-40 MeV protons for the different regions of the polar 
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caps. The configuration of these access windows will then be related to 
the magnetospheric models mentioned above, and hence to the question of 
magnetic field merging. 
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II. INSTRUMENT 

The detector system used in this study consists of two independent 
particle telescopes, each designed (a) to optimize charge and energy 
resolution within a range of incident particle energies, and (b) to be 
sensitive and operable over a wide range of incident particle fluxes. 
Because of their orientations on the spacecraft (see Section III), these 
telescopes are referred to as the vertical telescope and the horizontal 
telescope. Measurements consist of the counting rates, energy loss, and 
range of incident particles [28]. 

i 

The rapid motion of a polar orbiting satellite with respect to the 
geomagnetic field subjects a detector system on board to rapid fluctua- 
tions in incident fluxes, due to changing geomagnetic cutoffs and trapped 
particles. This tends to place an upper limit on the sampling period for 
the various rates monitored by the system. In order to obtain adequate 
resolution of these fluctuations, the rates for this experiment were 
averaged over a maximum of 288 msec, during which time the spacecraft 
will have moved, typically, 1.1 minutes in latitude. Although this 
averaging rate represents the situation for the bulk of the data re- 
ported here (data which were collected on the on-board tape recorder), 
when the spacecraft telemetry was in the real time mode the averaging 
time is either 72 msec or 18 msec, depending on the exact telemetry con- 
figuration. Data collected in the real time mode are particularly 
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suitable to the study of phenomena closely related to the geomagnetic 
field, such as cutoffs. 


Vertical Telescope 

Figure I I -1 shows a schematic cross sectional view of the vertical 
telescope, which consists of a stack of circular detectors and absorbers 
From top to bottom, they are: a gold-silicon surface barrier detector, 
an aluminum absorber, another gold-silicon surface barrier detec- 
tor, and a copper absorber. Such solid state detectors measure the 
energy lost by a charged particle passing through the sensitive region 
of the detector; the physical characteristics of these detectors and 
absorbers are given in table II-l. The detector stack is completely 
surrounded, except for the entrance aperture, by a cylindrical plastic 
scintillator cup, which serves to (1) collimate the response of the 
telescope, (2) help discriminate against side showers and nuclear inter- 
actions within the stack, and (3) determine the high energy analysis 
limit (^40 MeV/nucleon) for the system. The entrance aperture is 
covered by a 3/4 mil (0.00075 in) aluminized mylar light baffle; the 
thickness of this baffle contributes significantly to the low energy 
threshold (1.21 MeV/nucleon) of the telescope. 

Incident particles are analyzed only if not registered in the anti- 
coincidence detector D 3 . The number of particles with energy losses 
above an electronically determined threshold in Dj and D 2 are recorded, 
along with the number of DjD 2 coincidences. In addition, the energy 
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Figure II-l 


Schematic cross section of vertical telescope. 
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D, Au-Si DETECTOR 
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VALUES FOR OGO-EZ ARE IN PARENTHESES. 

D, AND D 2 BOTH HAVE DEPLETION 
DEPTHS OF 56 mg/cm 2 
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Table II-l 

Physical Properties of Detectors and Absorbers 
in the Vertical and Horizontal Telescopes 


Detector/Absorber 

Material 

Thickness 

(mg/cm 2 ) 

Sensitive 

Area 

(cm 2 ) 

Window (Vertical) 

Mylar 

2.6 ± 0.2 

- 

D (Total) 

1 (Depletion region) 

Silicon 

106 ± 2 
56 ± 5 

2.4 (nom.) 

Absorber 

Aluminum 

36 ± 4 

- 

n (Total ) 

2 (Depletion region) 

Silicon 

133 (nom.) 
56 (est.) 

3.5 (nom.) 

Absorber 

Copper 

1790 ± 45 

- 

Inner housing 

Magnesium 

138 ± 9 

- 

d 3 

Plastic 505 ±27 

Scintillator 

- 

Window (Horizontal) 

Mylar 

1.22 (nom.) 


u (Total) 

1 (Depletion region) 

Silicon 

11.7 (nom.) 
5.8 (nom.) 

0.079 (nom.' 
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lost in Dj is measured for one event during each sample period. The 
number of particles measured by D 3 are also recorded. The electronic 
configuration and the geometrical factors for the various rates are dis- 
cussed in more detail below. 

As noted instable II-l, the detectors used in this experiment are 
not fully depleted; fully depleted detectors with such large sensitive 
areas were not available when this experiment was designed. As will be 
noted below, the non-uniformity of the depletion depth over the sensi- 
tive area of the detector was one of the principle factors which limited 
the resolution of the pulse height analysis of the output from Di. 

Eorizontal Telescope 

A schematic cross section of the horizontal telescope is shown in 
figure II-2. This telescope, which is considerably smaller than the 
vertical telescope, consists of a single gold-silicon surface barrier 
detector completely surrounded, except for the entrance aperture, by a 
magnesium shield. As in the vertical telescope, the entrance aperture 
is covered by a 0.00035 in. aluminumized mylar light baffle. The physi- 
cal characteristics of the horizontal detector and window are included 
in table II-l. It should be noted that the collimation for this tele- 
scope is passive, in contrast to the active collimator available for the 
vertical telescope. The effect of this difference on the response of 
these two telescopes is indicated below in connection with the energy 
dependent geometrical factors. All charged particles losing more energy 
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Figure I I -2 


Schematic cross section of horizontal telescope. Note that the scale 
used here is different than that used in figure II-l. 
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in the sensitive region of the detector than the electronically deter- 
mined threshold are counted. 

Electronics and Data Formatting 

The block diagram in figure I I -3 illustrates the logic configuration 
of the electronics system associated with the two charged particle tele- 
scopes and the format in which the data are stored in the spacecraft. 

The details of the electronic components are similar to those for Exper- 
iment F-20 on 0G0-6 which were specified in a paper presented to the 
Fourteenth Nuclear Science Symposium of the IEEE in 1967 [29]. As indi- 
cated on this figure, the following information is stored by this exper- 
iment: 

1. The 2°, 2 4 , and 2 9 bits of the recycling scalars corresponding 
to the number of counts from the following four logic configu- 
rations: VjVT. V 2 V 7, ViV z V7, and ^ (VisDi,*V 2 sD 2i V 3 hD 3 ). 

2. The 2 4 , 2 5 , and 2 6 bits of a recycling scalar associated with 
the number of V 3 counts. 

3. A 256-channel pulse height analysis of the amount of energy 
deposited in D x . Only one such analysis can be stored; the 
contents of word 17 are erased prior to storing a new event 
analysis. Thus only the last event analyzed before each read- 
out is available. 

4. Flags indicating whether the pulse height analyzed event is a 
new event (since the last telemetry readout) and whether the 
threshold of 0 2 was exceeded on the event analyzed. 



Figure I I -3 


Functional block diagram of electronics system. 
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5. A digitization of an analog rate meter connected to the output 
of the V 3 discriminator. 

The spacecraft telemetry will be discussed in Section III. 

Thresholds and Calibrations 

The electronic thresholds associated with the discriminators of the 
four detectors were set so as to minimize the contamination due to back- 
ground noise. The values of these thresholds in terms of detector out- 
put and of incident particle energy were determined by electronic and 
particle calibrations. The electronic calibrations give a precise deter 
mination of the discriminator thresholds and pulse height analyzer re- 
sponse as a function of charge at the input to the charge sensitive pre- 
amplifier, while the particle calibrations enable a determination of the 
response and resolution of the detectors. 

In order to interpret these calibrations in terms of the response 
of the telescope, interpolations were made in the range-energy loss 
tables given by Janni [30] wherever necessary. These tables were gen- 
erated by integrating a semi -empirical expression for the energy loss of 
a charged particle passing through a homogeneous material. The cali- 
brations are described in more detail in a Space Radiation Laboratory 
Internal Report [31]. 

The results of the electronic calibrations of the detector thres- 
holds are given in table II-2. 



18 


Table II-2 

Energy loss thresholds for 0G0-4 detectors 


Detector 

Threshold 
(25 °C) 

Temperature 

Coefficient 

Di 

424±22 keV 

0.0 keV/°C 

d 2 

253±23 keV 

0.3 keV/°C 

d 3 

53±3 mV 

-0.1 mV/°C 

Hi 

377±14 keV 

0.6 keV/°C 



/, 
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By using the range-energy loss tables mentioned above, the corre- 
spondence between the pulse height analyzer channel thresholds and inci- 
dent proton energy can be determined. This is shown in figure I 1-4. An 
important facet of the correspondence shown in this figure is that for 
channel numbers greater than 29 (and less than 72) two incident proton 
energies can be associated with each channel threshold. Below channel 
29, the availability of the range information provided by D 2 makes the 
association between incident proton energy and channel threshold unique. 
The effect of these double-valued, pulse height analyzer channels on the 
analysis of the data will be discussed in Section IV. 

The proton response of the vertical telescope was determined by ex- 
posing the assembled telescope to a series of proton beams using the 
Caltech Tandem Van de Graff accelerator [31]. 

The electron response of the detectors in the vertical telescope 
was determined by exposing nearly identical detectors to the beam of a 
magnetic spectrometer [32]. Electron geometrical factors as a function 
of energy for Di and D 2 determined in this manner are given in figure 
I I -5 . It is clear from these results that the electron sensitivity of 
the ViV 2 V7 rate is negligible. The thin depletion depth and high elec- 
tronic threshold of the H x detector insure that the electron sensitivity 
of the H x rate is < 10‘ 4 [33]. 

The incident energy ranges corresponding to the rates measured by 
this experiment are summarized in table 1 1 -3 . In addition to the aver- 
age geometrical factors included in table I 1-3, energy dependent proton 
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Figure I I -4 

Correspondence between pulse height analyzer channel thresholds and inci- 
dent proton energy (incident in mylar window). The D x discriminator 
level is indicated, along with the energies corresponding to the trigger- 
ing of D 2 and D 3 . The alpha particle response is also shown, Note the 
set of channels for which there is no one-to-one correspondence between 
channel threshold and incident proton energy. These channels are referred 


to as double-valued. 
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21 



INCIDENT ENERGY(MeV/nucleon) 
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Figure I I -5 


Electron geometrical factors for VjV 3 and V 2 V 3 rates. 
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Table I I -3 

Incident Energy Ranges Corresponding to 
Electronic Rate Configurations 


Rate 

Configuration 

Incident Energy Range 

Electrons Nuclei 

(MeV) (MeV/nucleon) 

Approximate 

Geometrical Factors (Aft) 
(cm 2 -sr) 

Electrons Nuclei 

ViV 3 

0.45-1.8 

1.24-40.4 

0.08 

1.06 

v 2 vl 

0.67-2.5 

9.3-40.4 

0.36 

1.42 

ViV 2 VT 

0.67-1.8 

9.3-40.4 

10" 3 

1.16 

Hi 

— 

0. 88-^4. 5 

10' 4 

0.013 

v 3 

0.53 

37. 

? 

? 
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geometrical factors have been calculated for all rates but V3 using a 
Monte Carlo technique. These are shown in figures 1 1 -6 and 1 1 -7. Com- 
paring these two figures, the effect of the passive collimation in the 
horizontal telescope is clear: the high energy response of the H 2 rate 
is quite complex. 



Figure I 1-6 


Energy dependent geometrical factors for the V 2 V 3 , and V 3 V 2 V 3 rates. 
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PROTON ENERGY INCIDENT ON INSTRUMENT (MeV) 



Figure I I -7 


Energy dependent geometrical factor for the H! rate. 
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III. SATELLITE 

0G0-4 is the fourth in the series of six Orbiting Geophysical 
Observatory satellites sponsored by the National Aeronautics and Space 
Administration. This satellite was launched from the Pacific Missle 
Range on 28 July 1967 into an orbit with an initial inclination of 86°, 
apogee of 908 kilometers, perigee of 412 kilometers, and orbital period 
of 98 minutes [34]. 

The instrument described in Section II was mounted on the space- 
craft so that only particles incident on the vertical telescope from 
directions within 30° of the zenith do not trigger the anti -coincidence 
counter, D 3 , prior to passing through the solid state detectors, D x 
and D 2 (see figure II-l). Similarly, the horizontal telescope was 
mounted so that collimated particles represent particles incident from 
directions nearly perpendicular to the zenith. 

Because of a need to time-share the spacecraft telemetry facilities 
between two sets of experiments, the data collected by this experiment 
were only available about 50% of the time: often in two day periods 
separated by two day gaps in the data. The effect of this telemetry con- 
figuration on the data from this experiment is illustrated in the 0G0-4 
Data Coverage Plots [35]; figure III-l is a typical example of these 
plots . 
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Figure III-l 

A typical example of the 0G0-4 Data Coverage Plots. An explanation 
of the information displayed on this figure is given by Evans [35], 
where a complete set of these plots covering the entire period 
during which at least one of the 0G0-4 tape recorders was operable 
will also be found. These plots are designed to fulfill the dual 
purpose of indicating the availability of the data from this experi- 
ment and of comparing these data to other geophysical, interplanetary, 
and solar data of interest. 
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The spacecraft was deactivated on 23 October 1969, but the failure 
of the second of the two on-board tape recorders on 19 January 1969 
represented, at least for this study, the practical limit of the period 
for which useful data are available. 

The trajectories of 0G0-4 during several typical orbits mapped 
onto a polar representation of geocentric latitude vs. geocentric 
longitude are shown in figure I I 1-2 . As the plane of the satellite 
orbit remains relatively stationary in space, the earth rotates under 
the satellite, causing the apparent shift in the trajectories. Contours 
of constant invariant latitude (a) [36] are also projected into this co- 
ordinate system for reference; it should be noted that some of the orbits 
reach a much lower maximum invariant latitude than others, especially in 
the south. 

Since many of the charged particle effects which can be measured 
by this experiment are closely related to the geomagnetic field and the 
sun-earth-satellite orientation, the data observed with 0G0-4 can be 
more efficiently organized in a coordinate system reflecting the position 
of the satellite in the geomagnetic field and the sun-earth-satellite 
orientation. Of several coordinate systems reflecting these parameters, 
one of the most effective and widely used is that consisting of in- 
variant geomagnetic latitude (a) vs. magnetic local time (MLT) [37]. The 
trajectories shown in figure III-2, mapped into this system, are shown in 
figure I I 1-3 . It is clear that in this coordinate system the effects of 
the rotation of the earth are minimized, , and the differences in maximum 



Figure I I 1-2 


A projection of typical 0G0-4 trajectories onto a polar representation 
of geocentric latitude vs. geocentric longitude. Projections of these 
trajectories onto both geographic poles are shown for comparison, and 
contours of constant invariant geomagnetic latitude (a) are indicated. 
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Figure 1 1 1 -3 

A projection of the same 0G0-4 trajectories shown in figure I I 1-2 into 
a geomagnetic coordinate system: invariant latitude (a) vs. magnetic 
local time (MLT). The field coefficients given by Cain,et at. [38] 
were used to calculate a. 
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invariant latitude are depicted more explicitly. 

Since for the bulk of the observations reported here a critical 
distinction will be made between data collected at low and high polar 
geomagnetic latitudes (see section V), the data collected during polar 
passes such as 2039 and 2040 in the north and 2038, 2039 and 2040 in 
the south may be seriously degraded due to the limited degree to which 
the satellite penetrates to high geomagnetic latitudes. 
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IV. DATA ANALYSIS 


The processing scheme used for the data from this experiment is 
illustrated by the block diagram shown in figure IV-1. The data came 
to Caltech stored on two classes of magnetic tapes: experimenter tapes 
and attitude-orbit tapes. The former contained the decommutated data 
for this experiment, while the latter provided the position and orien- 
tation of the spacecraft as a function of time (one point per minute) 
calculated from data collected by tracking stations. The information 
from these two classes of tapes was interleaved, time ordered and 
stored on a third class of magnetic tape, referred to as an abstract 
tape. The study being reported here has involved the processing of 
data from a total of 931 magnetic tapes, representing almost 300,000 
minutes 62 * 10 6 read-outs) of playback data. The format of the^ 
data on these tapes is described elsewhere [39]. 

A complete collection of plots showing the detector rates as a 
function of time [28] serve as a catalogue of these data from which one 
can evaluate the performance of the experiment and specify the periods 
of data most suitable for a particular study. With this information any 
of the other programs indicated on figure IV-1 can be employed. 

The calculation of rates for the output of the solid state de- 
tectors is complicated somewhat by the fact that only the 2°, 2 4 and 2 9 
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Figure IV-1 


Block diagram representation of the procedure developed to process 0G0-4 
data at Caltech. The symbolism used is as follows: 



= magnetic tape 



= computer program (program names are given in capital 

letters.) 



= branch (decision) 



= cards punched 


= results (only the resultant plots are shown on this 
diagram; each program has the option of 
tabulating results as well.) 


The rate plots, which constitute the data catalogue, are discussed by 
Evans, et at. [28]. 
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The decision at this point depends en the type of 
problem being studied and the available data as Indi- 
cated on the rate plots and 0G0-4 Data Coverage Plots. 

Some of the choices currently available are Indicated. 
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bits of the rate scalers are recorded. Figure IV-2 illustrates that 
special care must be exercised in the calculation of rates in some 
ranges, and that for detection rates above ^500 counts per read-out 
(^1740 counts per second for playback data) an unambiguous determina- 
tion of these rates can be made only in special circumstances. 

o 

Using the accumulation of many pulse height analyzed events along 
with the response curve shown in figure I 1-4 and the energy dependent 
geometrical factor illustrated in figure II-6, differential flux spectra 
can be calculated. However, as pointed out in connection with the D x 
response curve (figure 1 1-4) , for a range of pulse height analyzer chan- 
nels the correspondence between channel threshold and incident proton 
energy is not one-to-one. As a consequence, the spectra presented in 
Section V have been calculated using the iterative technique illustrated 
in figure IV-3. 



Figure IV-2 


Probabilities of a change occurring in the 2 °, 2 4 and 2 9 -bits of a 
detector scaler during one read-out as a function of the event rate 
expressed in counts/read-out. 
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Figure IV-3 

Iterative technique used to calculate differential energy spectra. 

A spectrum is first calculated in a straightforward manner, and standard 
regression techniques are used to determine the "best" fit to 

dJ = AE _y 

using only the points for which the energy identification is unique 
(figure IV-3a). This function is then used to determine the proportion 
of events in a double valued energy bin to be assigned to each energy 
interval corresponding to the group of pulse height analyzer channels 
involved. The flux for this group of channels is then recalculated. 

The results of this procedure (figure IV-3b) are valid (for the affected 
points) only to the extent that the actual spectrum agrees with the 
function used to approximate it. 
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V. OBSERVATIONS 

Utilizing the data processing procedure outlined in Section IV, 
many of the available data from 0G0-4 have been processed. The 0G0-4 
Data Coverage Plots [35] indicate the periods during which data were 
available, and also indicate for which of these periods the data have 
been processed. These data are particularly appropriate to the study 
of the entry of solar particles into the earth's magnetosphere, and 
several periods have proven to be instructive. In addition to electron 
observations, two examples of proton observations will be presented as 
illustrative of the type of observations available. A more complete 
compilation of pertinent data, including those events presented here, 
is given in Appendix A. 

Electron Observations 

Although the instrument described in Section II was not designed 
primarily to detect electrons, several periods of high electron fluxes 
were observed in the polar cap regions on board 0G0-4. Some of the 
characteristics of electron polar cap observations have previously been 
reported [27,40,41], but the data presented here provide, for the first 
time, a comprehensive mapping of the boundary of the electron polar cap 
region. As discussed in Section VII, this boundary represents the 
boundary between open and closed geomagnetic field lines. 
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As mentioned in Section II, the VjV^ and V 2 VJ rates are much less 
sensitive to electrons than to protons. Although the contribution of 
electrons to the observed rates varies considerably, it is typically * 15 % 
for the ViTJ rate and *35% for the V 2 V 3 rate. In addition, we recall 
that the 0G0-4 anti -coincidence scintillator, V 3 , is sensitive to elec- 
trons above a threshold energy of ^530 keV. Normally, then, electrons 
constitute a rather minor constituent of the observed rates. There are 
periods, however, when the electron flux is sufficiently high relative 
to the proton flux that electrons become an identifiable constituent of 
the V 2 V 3 and V 3 rates. The identification of these periods and 

the identification of the rates as predominantly due to electrons is il- 
lustrated in figure V-l. The rates from this polar pass show a uniform 
electron polar cap flux between 1750:40 UT and 1805:40 UT; beyond 1806:00 
the spacecraft had moved to invariant latitudes where the electrons did 
not have free access. That this enhancement can be associated strictly 
with electrons can be seen clearly by comparing the V 2 VJ rate (electrons 
and protons) with the ViV 2 V 3 rate (protons only). In order to improve the 
precision with which the electron polar cap boundary could be specified, 
the V 3 rate was plotted on an expanded scale, shown in figure V- 2a. * 

Figure V-2b illustrates that the edge of the electron polar cap 
was normally associated with an electron spike near magnetic local 
midnight (MLT = 2100-0211). Figure V- 2c is an example of an electron 
polar cap observation showing a particularly striking feature (i.e. the 
sharp depression at 0233:40 UT). These features were rarely observed 
and normally occurred only late during the recovery phase of a magnetic 
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Figure V-l 

Rates observed during a typical 0G0-4 polar pass (a* 50°) illustrating 
an enhanced electron polar cap flux. These data are from the north, 
polar pass of Rev. 6291 on 26 September 1968. 
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Figure V-2 

High resolution plots of the V 3 rate, showing characteristic electron 
polar cap observations: 

V-2a: V 3 rate shown in figure V-l. 

V-2b: An example of the type of electron spike associated with the 
edge of the electron polar cap between MLT ^ 2100 and 
MLT ^ 0200. 

V-2c: An example of a well-defined feature in the electron polar 
cap flux. Such features were rarely observed. 
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storm. 

A total of 333 such observations were made during magnetically quiet 
periods, and the resultant map of the electron polar cap boundary is 
shown in figure V-3. Here the geomagnetic coordinates (invariant lat- 
itude and magnetic local time) of each observation of the boundary are 
indicated by a symbol. The apparent dependence on geomagnetic activity 
will not be discussed here. The values reported by Vampola [27], which 
he extrapolated to magnetically quiet conditions, are shown for compar- 
ison. 

We can now define the following terms with respect to the electron 
results from this experiment: high polar latitudes (PPL) and low polar 
latitudes (LPL) . High polar latitudes will refer to the invariant 
latitudes between the electron polar cap boundary and the geomagnetic 
pole; low polar latitudes will refer to the invariant latitudes be- 
tween geomagnetic cutoff for 1.2-40 MeV protons (below which polar 
proton fluxes cannot be observed due to the Earth's magnetic field) 
and the electron polar cap boundary. These definitions will be useful 
in organizing the proton observations. 

Proton Observations 

Most of the data collected by this instrument are by design 
dominated by proton fluxes. Since we are interested in studying those 
periods during which the proton counting rate changes significantly 
during a time scale of several hours, the proton data can be effectively 
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Figure V-3 


Observed electron polar cap boundary data in an invariant latitude 
magnetic local time coordinate system. 
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displayed on a plot of proton counting rate (or flux) as a function of 
time, which is referred to as a profile. In particular, we will be in- 
terested in comparing the profiles of fluxes measured at low polar lat- 

I 

itudes and at high polar latitudes. As discussed in Section VII, the 
north and south low polar latitude regions are most likely connected by 
closed field lines, while the north and south high polar latitudes are 
not. Therefore, no distinction will be made between data collected in 
the north and south low polar latitude regions, whereas the data from 
the two high polar latitude regions will be kept separate. 

The temporal variations observed in the proton flux can normally 
be divided into two main types: those events associated with and having 

the characteristics of solar flare events, and those events which have 
been variously referred to as Energetic Storm Particles (Bryant, et al. 
[42,43] and Rao, et dl. [44]), Delayed Particle Events (Lin and Ander- 
son [45]), and Protons Associated with Centers of Solar Activity (Fan, 
et al. [46]). We will follow the lead of Anderson [47] and try to avoid 
any semantics problem by referring to these events as edp (energetic 
Delayed Particles) events. Although treated in more detail in Section VI, 
a brief description of the essential differences between these two classes 
of events may indicate why the distinction between them is made. Flare 
events are characterized by impulsive ejection of particles by the sun 
followed by their propagation through interplanetary space. This typi- 
cally results in a rapid rise (^few hours) to a maximum proton flux, 
followed by a long, nearly exponential decay with a time constant of about 
one day. EDP events, on the other hand, are characterized by a region of 
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limited radial extent and enhanced low energy particle flux being 
convected outward by the solar wind plasma (expanding solar corona). 

For our purposes here, EDP events are phenomenologically characterized 
by the following (after Anderson [47]): 

1. Fluxes occur predominately at low energies. 

2. Profiles show much more rapid rise and fall than flare event 
profiles and often exhibit large fluctuations during the event. 

3. The duration of the event is normally between o(l hour) and 
°(1 day). 

4. Events are often associated with a "weak depression" in the 
sea level neutron monitor rates. 

5. Events are sometimes associated with one or more geomagnetic 
sudden commencements or sudden impulses. 

Presented below are typical examples of the observations of both 
of these classes of particle events. 

EDP event of 1 December 1967 

i 

Figure V-4 shows the VjVJ rate as a function of time for each 
available polar pass from 1700 UT on 1 December 1967 to 0700 UT on 
2 December. This figure illustrates the averaging intervals used to 
obtain average counting rates for each of the polar cap regions (LPL 
and HPL) , and the relation between these averaging intervals and the 
electron polar cap boundary data presented in figure V-3. The profiles 
of the flux measured in the low polar latitude (LPL) region and in the 



Figure V-4 


V X V7 counting rate (^15 second averages) during all of the available 
polar passes from 1700 on 1 December 1967 to 0700 on 2 December. The 
averaging intervals used to obtain average rates for low polar latitude 
( ■■■■■■■ ) and high polar latitudes ( r~ ~ i ) are shown and 

compared to the approximate location of the electron polar cap boundary 
(a) as shown in figure V-3. In addition to universal time, geomagnetic 
coordinates (a and MLT) are also indicated. 
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north and south high polar latitude (HPL) regions are in figure V-5. 

All data points shown in this figure have errors comparable to those of 
the four points with explicit error bars. The arrows near the top of 
the figure indicate a gap in the data, which included both a south and 
a north polar pass (see Section III). The time resolution of the LPL 
profile is twice that of the other two profiles because the north and 
south LPL regions are considered equivalent. 

The identification of this event as an EDP event is supported by 
the extremely fast "decay" of the flux, the predominantly low energy 
nature of the event (see discussion of spectra below), and the presence 
of a weak depression in the Alert Neutron Monitor [48]. In addition, 
figure V-6 shows a direct comparison between the data in figure V-5 and 
the interplanetary flux of 0.79-9.6 MeV protons and 0.17-1.00 MeV 
electrons measured by the University of Chicago experiment on board IMP-F 
[49]. The correlation between these interplanetary data and the LPL pro- 
file is notable. 

The most obvious features of these profiles are that (1) all three 
profiles are different (a peak occurs first at LPL, next at northern HPL, 
and last at southern HPL) and (2) only the LPL peak occurs at the same 
time as the peak seen in the IMP-F data. 

In order to determine the relationship among these peaks, spectra 
were calculated at eight times during this event; the data for which 
spectra were calculated are indicated by the circled points on figure V-5. 
These spectra are shown in figure V-7. The data upon which the spectra 
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Figure V-5 

Flux profiles for EDP event of 1 December 1967. The profile for each 
region is distinguished from the other profiles according to the 
following code: 

• • • low polar latitudes 

• north high polar latitudes 

• • south high polar latitudes 



N-HPL 


\ - 

\ © 
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♦ S-HPL 

DATA POINT AT 
WHICH SPECTRUM 
WAS CALCULATED 
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Figure V-6 

A direct comparison between interplanetary and polar cap particle obser- 
vations for the EDP event of 1 December 1967. Interplanetary data are 
the flux of 0.79-9.6 MeV protons and 0.17-1.00 MeV electrons measured 
by the University of Chicago experiment on board IMP-F [49]. Polar cap 
data are the flux of 1.2-40 MeV protons and 0. 4-1.0 MeV electrons meas- 
ured by VjVJ on board 0G0-4. The interplanetary data have been multi- 
plied by 0.32 to normalize them to the 0G0-4 flux at the EDP peak. 
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Figure V-7 

Differential energy spectra of the following fluxes calculated from 
data collected during the 1 December 1967 EDP event: 

V-7a: Low polar latitudes prior to the beginning of the event. 

V-7b: Low polar latitude peak, northern high polar latitude peak, 

and southern high polar latitude data from the same time. 
V-7c: Flux from each of the three regions collected at or near 

the southern high polar latitude peak. 

V-7d: Low polar latitudes after the end of the event. 

In all four cases, LPL data points are indicated by an open circle (o), 
northern HPL by a closed circle (•) and southern HPL by a cross (x). 
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in figure V-7b were based were collected during consecutive polar passes 
in order to minimize any differences due to temporal effects; the same 
is true of the data in figure V-7c. We note the striking interchange in 
roles from figure V-7b to figure V-7c between the southern HPL region 
and the LPL and northern HPL regions. 

These eight spectra are summarized in table V-l , which gives the 
coefficient and exponent of the best fit for each spectrum to an equation 
of the form 

3F = ae " Y 

where J is the flux in parti cl es/cm 2 -sec-sr, and y is referred to as the 
spectral index. One standard deviation errors for A and y, as well as 
the x 2 "goodness of fit" parameter, are also indicated on this table. 

The spectra observed at each of the peaks indicate the following: 

(a) compared to the spectra observed before and after the event, each of 
the peaks consisted of a larger proportion of low energy protons and ex- 
hibited spectral indices typical of those observed for EDP events in inter- 
planetary space [50], and (b) the spectral index of the flux was essenti- 
ally the same at each peak. The first point supports the identification 
of these peaks as due to EDP fluxes, while point (b) is consistent with 
the interpretation that all of the peaks represent a sampling of the 
same particle population. This interpretation is aided by the obser- 
vations that the magnitude and spectrum of the ambient flux change very 
little throughout the event and that the spectral indices for the peak 
fluxes are significantly different than the index for the ambient flux. 



TABLE V-l 


1 Dec. 1967 EDP Event Spectra (see figure V-7) 
Best Fit to 4 f=AE _y 


Relationship to Event Coefficient 

"A" 

(cm 2 -sec-sr-MeV)“ 1 

Spectral 

Index 

n Y ii 

X 2 

P(>X 2 ) 

Pre-event LPL 

1.07 ± 0.18 

2.31 ± 

0.23 

9.33 

0.316 

LPL peak 

9.90 ± 0.15 

3.52 ± 

0.22 

0.96 

0.995 

Northern HPL peak 

10.90 ± 0.11 

3.67 ± 

0.17 

2.52 

0.961 

Southern HPL at LPL peak 

1.21 ± 0.21 

2.53 ± 

0.29 

1.47 

0.983 

LPL at Southern HPL peak 

1.13 ± 0.19 

2.66 ± 

0.24 

2.19 

0.974 

North. HPL at S. HPL peak 

1.10 ± 0.23 

2.52 ± 

0.32 

1.95 

0.962 

Southern HPL peak 

4.43 ± 0.26 

3.62 ± 

0.41 

1.28 

0.989, 

Post-event LPL 

1.47 ± 0.21 

2.69 ± 

0.29 

1.08 

0.993 
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These observations lead to the following description of the 1 Decem- 
ber 1967 event. The interplanetary region of enhanced flux arrives at 
the earth and is observed almost simultaneously by IMP-F and at LPL. 

About 1.2 hours later, by which time the LPL flux has decreased signif- 
icantly, this flux is observed to reach a maximum at northern HPL. 
Finally, after these two fluxes have definitely returned to ambient flux 
levels, the EDP flux is observed at southern HPL. 

Solar Flare Event of 2 November 1967 

The second of the two classes of solar particle events of interest 
here consists of those events associated with solar flares. As with EDP 
events, we will compare the profiles of the flux at low polar latitudes 
and the flux at north and south high polar latitudes. Many flare events 
have been used in this study; most of these are subject to the limitation 
imposed by the time-sharing nature of the 060-4 telemetry: in general, 
data are available for two day periods separated by two day gaps. This 
limitation is more critical for flare events, which might last 6-8 days, 
than for EDP events, which have typical durations between one hour and 
one day. 

Typical of the observations which, in spite of this limitation, 
have features of interest are the profiles of the 2 November 1967 flare 
event in figure V-8. Preliminary results of the study of this event have 
been previously reported [11]. Typical errors are indicated for arbitra- 
rily selected data points. 
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Figure V-8 

Profiles of the 2 November 1967 solar flare event. The symbols in this 
figure are explained in Appendix A (see table A-3). A 2B flare accurred 
on the sun at 0852 UT on 2 November 1967, accompanied by 2-12 $ X-ray 
emission that peaked at 0858 UT. 


FLUX (cm 2 -sec-sr) _l 
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This is a particularly good example of a prolonged difference be- 
tween the flux observed in one high polar latitude region and the flux 
observed in the other two polar regions. The north high polar latitude 
flux is consistently lower than the other two fluxes for a period of ^25 
hours. This type of profile configuration, where one high polar lati- 
tude profile is markedly different than the other high polar latitude 
profile has been referred to as a North-South (N/S) asymmetry [11,25]. 

The appearance of such a N/S asymmetry during a single polar pass 
is illustrated by figure Y-9. This figure shows the V 1 Y 3 rate for the 
north and south polar passes of one orbit early during the 2 November 
1967 event. Southern polar passes from later in the event show even less 
structure. The ratio of the flux of 10-40 MeV protons to the flux of 
1.6-40 MeV protons is an indication that the same ratio of minimum to 
maximum flux occurred over the entire measured energy range from 1.2 
to 40 MeV. As is discussed below, pulse height analysis gives the same 
indication. Figure V-9 also contains the ratio of the vertically inci- 
dent to the horizontally incident flux, normalized so that unity cor- 
responds to an isotropic flux. There is no evidence of a large anisotropy 
in conjunction with the intensity variations, ruling out a pitch-angle 
dependent cut-off effect. 

The manner in which the feature in the north polar proton flux il- 
lustrated in figure V-9 is repeated consistently for a long period as 
shown in figure V -8 suggests the possibility of referring to these obser- 
vations as persistent features. A feature in the proton polar cap flux 
will be termed persistent if it is observed during two or more 



74 


Figure V-9 

Selected north and south polar passes observed early during the 2 No- 
vember 1967 solar flare event. Isotropy is indicated by a vertical/ 
horizontal ratio of unity. Any spectral change would be indicated by 
a change in the 10-40 MeV/1.6-40 MeV proton ratio, which has a statis- 
tical error of $1%. Invariant latitude is included for reference. 
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consecutive orbits (excluding orbits for which appropriate data are 
unavailable). 

The mapping of the regions of maximum and minimum intensity into the 
invariant latitude-magnetic local time (a-MLT) coordinate system is shown 
in figure V-10, along with the low latitude rigidity cutoff previously 
determined for 1.5 MeV protons during a magnetically quiet period in 
1961 [51]. 

Spectra were calculated throughout this event, and figure V-ll 
shows a typical example of the spectrum from the northern high polar 
latitude region compared to the spectrum from the low polar latitude 
region observed at the same time. The shapes of these spectra are re- 
markably similar and differ only by a constant multiplicative factor. 
Although there are some variations early in the event, as shown on figure 
V-12, no prolonged energy dependent effects can be observed. 

Comparison of Electron and Proton Observations 

During the 28 September 1968 event, a comparison between electron 
polar cap observations and the features which occur in the polar proton 
intensities was possible, and it is shown in figure V-13. During this 
event the persistent proton feature took the form of a depression. The 
observations of the electron polar cap boundary are mapped into an invari- 
ant latitude-magnetic local time coordinant system and represented by solid 
circles; the high latitude limits of the enhanced flux region associated 
with low polar latitudes appear as open circles. It is clear that the 
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Figure V-10 

A mapping of the north polar passes of 0G0-4 into invariant latitude- 
magnetic local time coordinates for the period of the 2 November 1967 
flare event. The locations of the regions of maximum and minimum 
proton fluxes are shown. Observed geomagnetic cutoffs are indicated 
and compared to cutoffs reported for a quiet geomagnetic field in 
September, 1961 [51]. The extent of data coverage is indicated where 
it is less than a complete polar pass. 
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Figure V-ll 

Typical proton differential energy spectra observed during the 2 Novem 
ber 1967 solar flare event. Open circles represent the spectrum ob- 
served at low polar latitudes, and crosses represent the spectrum 
observed at the same time at northern polar latitudes. 







81 


Figure V-12 

Ratio of low polar latitude flux to northern high polar latitude flux 
as a function of energy for the period of the 2 November 1967 solar 
flare event. Each symbol consistently represents the ratio for a 
specific range of energies. The longer horizontal lines indicate the 
±lo limits of the ratio of the low polar latitude ViV7 rate (1.2-40 MeV 
protons) to the northern high polar latitude rate. 




REV. NO. DURING 2 NOVEMBER 1967 EVENT 
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Figure V-13 

Comparison in invariant latitude-magnetic local time coordinates of 
observed electron polar cap boundaries with observed high latitude limits 
of regions with enhanced proton fluxes. During the event shown, these 
enhanced fluxes were observed at low polar latitudes. 
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latitudes of the proton limits are * the latitudes of the electron polar 
cap boundary. This comparison indicates that the region where electrons 
have rapid access is the region where the proton access is delayed. 
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VI. BACKGROUND 

The significance of the data presented in Section V can be most 
effectively investigated in the proper context. As will be pointed -out 
in Section VII, these data are particularly relevant to questions con- 
cerning the configuration of the distant geomagnetic field. Since 
charged particles observed in the polar regions propagate through inter- 
planetary space, a brief description of this environment and the behavior 
of energetic solar particles therein is appropriate. In order to estab- 
lish the context for the discussions to follow, the distorted configur- 
ation of the near-earth geomagnetic field will be presented, along with 
the significance and difficulties of studying the distant geomagnetic 
field. 

Interplanetary Environment 

The model which we will use here is the one which was first proposed 
by Parker [52]. In this model the solar corona plasma continues to. ex- 
pand radially outward from the sun to form a super-Alfvenic solar wind. 

The solar magnetic field is "frozen into" the solar wind plasma [53] and 
is thus convected outward from the sun. In the presence of a homogeneous, 
uniform solar wind flow, the combination of the radial solar wind flow 
and the rotation of the sun would cause the solar magnetic field to be 
pulled into an Archimedian spiral configuration, such that near the earth 
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the angle between the magnetic field and a radial vector would be about 
48°, measured in the ecliptic plane [54]. Neither the source of the mag- 
netic field nor the medium in which it is being convected is either homo- 
geneous or uniform, however. As a consequence, the configuration of the 
interplanetary magnetic field fluctuates considerably from this ideal po- 
sition (Cf. [55]), and the actual configuration is more like that depic- 
ted schematically in figure VI-1. An important parameter of the solar mag- 
netic field is the mean direction of the field, designated as positive (nega 
tive) if the field direction Is predominately away from (toward) the sun 
along the average spiral direction. Most of the time this mean direction is 
observed to be divided into sectors [56], also shown in figure VI-1. The 
history of this sector structure during the lifetime of the 060-4 data ac- 
quisition [57,58] is indicated on the 060-4 Data Coverage Plots (an example 
is shown in figure IV-1 ; the complete set is given by Evans [35]). 

As indicated in Section V, solar cosmic ray events can be divided 

into two different phenomenological classes: (a) the prompt event, or 

solar flare event, in which particles arrive within a few hours of the par- 

ent flare, and (b) delayed events, or EDP (Energetic Delayed Particle) 
events, in which particles arrive a day or more after the parent flare [59]. 
Since each of these classes of events has been observed with 060-4, it is 
appropriate to briefly discuss some of the phenomenological aspects of 
each. 


Solar flare particles are thought to be injected into interplanetary 
space impulsively at the sun more or less coincident with a relatively 



Figure VI-1 


Schematic representation of the configuration of the solar magneti 
field near the ecliptic plane. 
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large release of energy (e.g. optical. X-rays, synchrotron radiation 
from energetic electrons). After injection, the propagation of these 
particles from the sun to the point of observation is determined by the 
characteristics of interplanetary space. The results of this propagation 
are normally an increase in the interplanetary flux to a maximum over a 
time scale of about five to fifteen hours [59], followed by a decay which 
becomes more or less exponential after a few hours with a decay time con- 
stant on the order of a day. The radial gradient in the flare flux is 
normally negative at 1 AU early in the event, but may become positive 
during the decay phase. Such a positive gradient has been observed by 
^ = 2.0 x 10" 3 %/Rq [60]. Current models for the pro- 
pagation of solar flare protons in interplanetary space [61,33] give 
j - 1.5-11.0 x io " 3 %/R^ as typical for gradients near 1 AU for late 
times (t > 140 hours). Early during the decay phase, interplanetary gra- 
dients could be as high as 40-50 x 10" 3 %/R^, depending on the position 
of the flare on the sun [33]. 

Anisotropies in the interplanetary proton fluxes associated with 
prompt events are observed throughout most events [59,62]. The degree 
and direction of the anisotropy varies greatly, however, during each event 
and also from one event to another. During the early phase of a prompt 
event, particles are normally observed streaming out from the sun along 
the direction of the interplanetary magnetic field [63,64]. This field 
aligned anisotropy decays rapidly and is usually not dominant shortly 
after the maximum flux is observed (5 to 15 hours after onset) [59,63]. 
After this time the anisotropy direction is independent of the magnetic 


* 

O' Gallagher 
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field direction and is inferred to be parallel to the solar wind velocity 
for a period of a few days [59]. The magnitude of the anisotropy during 
this phase is energy dependent and is normally between 5% and 20% [63]. 
Late in the decay phase of the event (more than four days after onset) 
the anisotropy is directed more or less perpendicularly to the magnetic 
field direction and its magnitude is -5% [65]. The general duration and 
magnitude of interplanetary anisotropies during these three phases will 
be referred to in Section VII when the effect of these anisotropies on 
polar cap observations is discussed. 

Another aspect of interplanetary anisotropies which will be of sig- 
nificance in Section VII is that, in general, the magnitude of the maximum 
anisotropy is inversely dependent upon the solar longitudinal separation 
between the parent flare and the foot of the interplanetary line of force 
passing through the point of observation [59]. Although this is generally 
true, large (> 3:1) anisotropies have been observed at lower proton ener- 
gies during the initial phases of some flare events which are time cor- 
related to east limb optical flares [66]. 

The other class of solar particle events consists of the delayed 
events, or EDP events, which have been observed many times in interplan- 
etary space (e.g., as cited in Section V, [42-47,50]). Most of the 
phenomenological aspects of these events are mentioned in Section V, so 
we will confine our attention here to their source, with one exception: 
these events are usually associated with large field directed anisotropies 
in interplanetary space [44]. These EDP phenomena are normally associated 
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with the so-called "active" regions of the sun; the enhanced flow of 
plasma from such a region sets up an interplanetary configuration such as 
that illustrated in figure VI -2 [cf. 42,46,47], which co-rotates with 
the sun. The effect of simultaneous observations at radially separated 
points in space would be to observe an abrupt and short-lived enhancement 
at the most sunward observation point first and at the point furthest 
from the sun last. The importance of this class of events to our study 
is that such events are associated with spatially well-defined inter- 
planetary features which sweep past the earth and down the geomagnetic 
tail (see below) at a nearly constant velocity, that of the solar wind. 
Polar cap observations of EDP events will provide the basis for determ- 
ining approximately where the particles observed in a given region of 

! 

the polar cap gained access to the geomagnetic field. 

Geomagnetic Field and Magnetic Merging 

An important consequence of the existence of the solar wind is the 
resultant configuration of the earth's magnetic field. The result of this 
plasma-magnetic field interaction is the confinement of the geomagnetic 
field to a cavity, called the magnetosphere, whose size and shape are de- 
termined by a balancing of the various magnetic and plasma pressures in- 
volved [67-70], The current picture of the near-earth configuration of 
the magnetosphere is illustrated in figure VI-3 [cf. [71]). As a matter 
of definition, those geomagnetic field lines on which some charged part- 
icles are capable of completing at least one bounce period (i.e., mirror- 
ing in both hemispheres) will be termed closed field lines. Those on 
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Figure VI-2 

Schematic representation of an interplanetary region of enhanced low 
energy flux co-rotating with the sun [42,46,47]. As the region co- 
rotates past the earth, radially separated observation points would 
record different event "arrival" and "departure" times. The differ- 
ence between these times would depend only upon the solar wind velocity 
and the radial separation between the points. 
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OBSERVATION POINT RADIALLY 
SEPARATED FROM THE EARTH 
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Figure VI-3 


Schematic representation of the near-earth configuration of the 
earth's magnetosphere. 



Magnetopause 



Magnetopause 
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which such trajectories are impossible for any charged particles will be 
termed open field lines. In general, the open field lines constitute the 

j 

field in the distant geomagnetic tail (see, for instance, [72,73]). The 
configuration of the geomagnetic tail has been mapped rather completely 
at least as far behind the earth as about 80 R^ (somewhat beyond the 
orbit of the moon) [71,72,74-76]. There is some indication that at 
least a geomagnetic wake may extend as far as 500 [77,78] or 1000 R^ 

[79-81] behind the earth. These observations are, however, complicated 
by the motion of the tail [80], and the detailed structure at these dis- 
tances has not been measured directly. 

A concept which will be of importance later should be defined at 
this time: that of an access window for the entry of particles into the 
magnetosphere. For a given location in the interior of the magnetosphere, 
the access window for particles of a given rigidity will be defined as 
the set of all points on the "surface" of the magnetosphere where particles 

i 

of that rigidity can gain access to the interior and be subsequently ob- 
served at the given location. Since a surface per se may be well-defined 
neither in the distant tail nor in a region where the two fields are 

i 

directly connected, the use of the term "surface" is meant to imply (a) 
the interface between the geomagnetic field and the interplanetary magnetic 
field and/or (b) the extension of this interface through a region where 
the fields are directly connected, if such a region exists. This defin- 
ition is similar to, but more formal than, that given by Gall, et at. for 
what they refer to variously as "penetration regions" [82] and "windows" 
[83]. 
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In the absence of direct measurements of the distant geomagnetic 
tail (beyond 80 R^), one of the currently unresolved questions concerning 
the geomagnetic field is whether or not this field merges (i.e. reconnects) 
with the interplanetary magnetic field. The rdle of charged particle ob- 
servations in the resolution of this question is that they can act as 
"probes" of a magnetic field where direct measurements of the field are 
infeasible. Assuming either the presence or absence of magnetic field 
merging between the solar and terrestrial fields, it is possible to post- 
ulate what the resultant magnetic field configuration might be. An 
analysis of a configuration so obtained will then yield implications 
concerning the characteristics of charged particle access from inter- 
planetary space to the interior of the magnetosphere and the relationship 
between interplanetary particle fluxes and fluxes observed in the polar 
cap regions. Comparisons of polar cap observations with the predictions 
concerning these observations arising from such a magnetic field con- 
figuration will yield constraints to be incorporated in the model. 

The implications of merging between the solar and terrestrial magnetic 
fields are not, however, restricted to questions of particle access into 
the magnetosphere. It is generally thought, for instance, that magnetic 
field merging may play a role in the formation of solar flares [5], and 
might even be involved in galactic phenomena [6]. An investigation of these 
possibilities is hampered, though, by the lack of concrete evidence con- 
cerning the conditions under which merging could take place. 


The basic mechanism for the reconnection of magnetic fields is 
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provided by the finite conductivity of the plasma. In 1956, Sweet [84] 
pointed out that the diffusion rate for magnetic fields in the vicinity 
of a neutral point (at which $ vanishes) will be much faster than in 
other regions of a plasma. Since the magnetic forces vanish at such 
a point, the magnetic field will tend to form a sheet near which 3 
reverses direction on a relatively small length scale; such a sheet 
is termed a neutral sheet. A neutral point also tends to be a stag- 
nation point of the plasma flow, in which case the magnetic field is 
convected toward the neutral point, where the field lines are "broken" 
and "reconnected" as the direction of the plasma flow changes by 90°. 

Since the field vanishes at the neutral point, however, the process is 
more accurately thought of as one in which the field being carried to 
the neutral point disappears and reappears in a different configuration, 
as determined by the fluid flow. Figure VI-4 illustrates this process, 
which is referred to as Sweet's mechanism; its astrophysical implications 
have been investigated by Parker [85]. 

In order to provide a mechanism whereby higher merging rates than 
those associated with Sweet's mechanism could be attained, Petschek [86] 
has suggested that while diffusion may be presumed to dominate the magnetic 
merging in the immediate vicinity of the neutral point, the plasma which 
is convecting the field away from the neutral point (after merging) may 
be moving at the Alfven velocity, V ft . Petschek finds that the maximum in- 
coming plasma velocity which can be supported by such a scheme is given 
by the following recursive relation: 
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Figure VI-4 

Schematic depiction of a neutral point, at which the magnetic field 
is shown undergoing merging via Sweet's mechanism. The angle between 
the two fields near the neutral sheet has been exaggerated. 
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where a is the conductivity of the plasma at the neutral point, and L 
is the scale size of the diffusion region. Sonnerup [87] has recently 
refined this approach somewhat, and has reported the following relation- 
ship for determining the maximum plasma velocity at which the fields may 
merge: 

U m = v a (1+/2) (6 ' 2) 

By treating the basic field equations involved, Yeh and Axford [88] 
conclude that in general there is no maximum merging velocity, although 
they show that Sonnerup' s solution is a special case of their solution — 
the only non-singular case, according to Sonnerup. Considering the 
scarcity of observational information concerning magnetic merging, it is 
virtually impossible to determine which, if any, of these solutions is 
more nearly correct. Even establishing the absence of merging between 
the solar and terrestrial magnetic fields would be significant. 

As a matter of terminology, geomagnetic field configurations arising 
from the assumption of the absence of magnetic merging between the solar 
and terrestrial fields are referred to as closed magnetospheric models, 
while those configurations based on the presence of merging are termed 
open magnetospheric models. In the next Section, after a discussion of 
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the general implications of the electron and proton observations, ex- 
amples of each of these types of magnetospheric models will be presented 
and compared to the 0G0-4 results presented in Section V and Appendix A. 
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VII. DISCUSSION 

I 

The information in Section VI provides the basic framework into 
which the data presented in Section V can be placed. The 0G0-4 electron 
observations provide a comprehensive determination of the boundary be- 
tween low polar latitudes (closed field lines) and high polar latitudes 
(open field lines). After this interpretation of the electron data has 
been established, the proton observations can be used to investigate the 
characteristics of the "windows" where these charged particles can gain 
access to the geomagnetic field. The observations of EDP events will be 
of particular interest in this regard, especially with respect to charac- 
terizing the positions and extents of these access windows. These EDP 
data represent the first reported observations of clearly identifiable 
interplanetary features, which are propagating down the tail, appearing 
at different times in different regions of the polar cap. With the par- 
ticle access configuration established, the three major models of the 
distant geomagnetic field will be presented, noting the relationship of 
each to the fundamental question of magnetic merging. The pertinent pre- 
dictions and extrapolations of the field configuration resulting from 
each model will be compared with the 0G0-4 data with the result that 

I 

severe constraints can be placed on each model. 
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Electron Polar Cap 

The low magnetic rigidity and high velocity of the 0. 4-1.0 MeV 
electrons observed by this experiment make these electron observations 
of particular importance in defining the transition between the region 
of open field lines and the region of closed field lines in the magneto- 
sphere [7-9,27,41,89,90]. The identification of the electron polar 
cap with open geomagnetic field lines can be established as follows. 

The origin of 50 keV polar cap electrons has been shown to be beyond 

about 64 R behind the earth by moon shadowing experiments conducted with 
ffi 

a lunar orbiter by Van Allen and Ness [89], Van Allen [7,8], Anderson 
and Lin [9], and Anderson [90]. In addition. Van Allen and Ness [89] 
observe a diffusion rate perpendicular to the magnetic field less than 
100 km/sec between the moon and the earth. These observations, coupled 
with the low rigidity of the electrons, imply that the edge of the elec- 
tron polar cap plateau should reflect the boundary of open field lines 
closely. The high electron velocities mean that differences in mode or 
position of access which might cause large variations in proton fluxes 
are not significant in electron observations. For instance. Van Allen 
[8] reports transit times for electrons travelling from interplanetary 
space to inside the magnetotail of <100 seconds, which were interpreted 
by Van Allen as indicating that electrons can gain direct access to the 
geomagnetic tail. As a result of the high velocities and rapid access, 
the electron flux over the polar caps is relatively homogeneous above the 
electron polar cap boundary (see figures V-l and V-2, also Vampola [27] 
and West and Vampola [41]). The relative uniformity of the electron flux 
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inside the boundary makes the specification of the location of this 
boundary rather unambiguous. We therefore interpret the data presented 
in figure V-3 as a definition of the region of the polar cap which is 
associated with open field lines in the geomagnetic tail. Whereas Vam- 
pola [27] has reported 25 observations of this boundary which he has ex- 
trapolated to geomagnetically quiet conditions (indicated on figure V-3), 
the 333 boundary observations made with 0G0-4 during magnetically quiet 
periods constitute a much more comprehensive mapping of the open field 
line region than has previously been possible. 

For the purpose of clarifying the proton observations, it is clear 
that the region defined as the high polar latitude (HPL) region in 
Section V is associated with the geomagnetic field lines which constitute 
the geomagnetic tail. Likewise, the low polar latitude (LPL) region is 
associated with closed geomagnetic field lines, which never extend far 
from the earth (< 30 R^) . In order to avoid the complications associated 
with any transition region, data associated with high polar latitudes 
have been collected at invariant latitudes well above the boundary in- 
dicated on figure V-3. Low polar latitude data similarly avoid the region 
close to the boundary of open field lines. As shown by figure V-18, per- 
sistent proton features can be identified consistently with the region 
of open field lines. With this in mind, we turn now to the proton obser- 
vations. 
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Access of 1.2-40 MeV Protons 

While electrons are able to define the edge of the open field line 
region with some precision, the much lower velocities and higher rigid- 
ities of the 1.2-40 MeV protons make them more sensitive to the effects 
of interplanetary gradients and anisotropies, different access modes and 
access window positions. Consequently, they can be used as "probes" to 
investigate the characteristics of the access windows. 

Prior to the report of preliminary results from this experiment for 
the 2 November 1967 event [11], some evidence of polar cap features was 
available [12-20], but the observations were too sparse to lead to com- 
prehensive analysis. Although since the report of the persistent feature 
observed during the 2 November 1967 event several observations of persis- 
tent features have been reported [21-26], the conclusions drawn from these 
observations still suffer from a paucity of data: only Bostrom [21] and 
Morfill and Quenby [26] have reported observations from more than one 
event (two events in both cases). 

A summary of all of the proton events observed to date in the 060-4 
data is presented in figure VII-1. As indicated in this figure, there is 
a strong correlation between the pole in which a persistent feature was 
observed on 0G0-4 and the interplanetary sector structure. With notable 
exceptions, the persistent features were observed in the north pole during 
negative sectors and in the south pole during positive sectors. This 
correlation is consistent with other reported studies in which the inter- 
planetary sector structure has been considered [13,24-26,91-93 (the 
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Figure VII-1 

Summary of 0G0-4 observations of solar events. The information indicated 

on this figure is as follows: 

1. Each event is identified by date. 

2. The phase of the flare during which a persistent feature was observed 
is indicated: 

/ = observation during the rise of flare 
\ = observation during the decay of flare 
E = EDP event observation 

3. The pole in which the persistent feature was observed is indicated 

( “" ) and compared to the interplanetary magnetic field sector 

LTV 

4. The duration (in hours) of the feature is indicated. In 75% of these 
observations the actual duration of the event was not measurable due 
to incomplete data coverage. 

5. The history of the north/south component of the solar magnetic field 
is represented by a vertical sequence of Vs (northward component) 
and Vs (southward component). Each symbol represents the average 
sense of this component during one hour [94]; the lowest symbol 
indicates the status of this component at the beginning of the 
observations . 
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Figure VII-1 (continued) 

6. Those events during which an enhancement was observed at high polar 
latitudes are indicated (*). 

7. For each polar pass of the indicated pole during the event, the ratio 
of the LPL ViVj rate to the HPL VjVj rate was calculated. The maxi- 
mum of these ratios is plotted as the "maximum relative 'depth' of 
the feature." 
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last three do not involve satellite proton polar cap observations)]. In 
accordance with this correlation and for the sake of clarity, the following 
terms will now be introduced: a-pole and 3-pole. These terms have been 
carefully selected in an attempt to avoid any terms which might connote 
any feature of any model . The correspondence between the north and south 
geomagnetic poles and the a-pole and 3-pole is dependent upon the inter- 
planetary sector and is given in table VII— 1 • With this definition, the 
field lines in the geomagnetic tail associated with the a-pole are usually 
more nearly anti parallel to the interplanetary magnetic field lines than 
those field lines in the geomagnetic tail associated with the 3-pole. We 
can now restate the sector correlation noted above as: "With notable ex- 

ceptions, the persistent features were observed in the 3 -pole." The per- 
sistent flare event features were correlated with the interplanetary sec- 
tor structure in this manner 91% of the time, while they appeared to be 
completely uncorrelated with the north-south component of the interplane- 
tary field. The interplanetary field had a southward component only 54% 
of the time, and figure VII-1 shows that during several events this compo- 
nent changed direction several times with no notable effects. 

The data in table VII-2 tend to strengthen the correlation between 
the interplanetary sector structure and the observations of persistent 
features considerably. This table lists all of the instances when a sec- 
tor reversal was observed in interplanetary space [94] while a persistent 
polar cap feature was being observed. Within the constraints imposed by 
the inherent time resolution of the observations and the frequent gaps in 
the data, the following is implied by these data: persistent features 



112 


TABLE VII-1 


Correspondence Between a-pole/ e-pole 
North/South Geomagnetic Poles 



| North 

South 


Pole 

Pole 

Positi ve 
Interplanetary 

a-pole 

e-pole 

Sector 



Negative 

Interplanetary 

g-pole 

a-pole 

Sector 
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TABLE VI 1-2 

Time Correlation of Sector Reversals [94] with Persistent Features 


Date 

Last Feature Obs. 

Hours 1V3 

Before (cts/sec) 
SR LPL HPL 

Sector 
Reversal 
MJni v . 
Time 
(HHMM) 

1st Pass after 
SR in Same HPL 
Hours 'W jV7 
after (cts/ 

SR sec) 

Figure 

Showing 

Profile 

11 Aug 67 

0.2 

28. 

40. 

1600 

1.5 

26. 


28 Oct 67 

0.6 

0.8 

0.4 

1200 

1.0 

0.55 


3 Nov 67 

1.3 

89. 

28. 

1300 

1.8 

45. 

E-3 

10 Nov 67 

0.9 

3.4 

2.2 

2300 

4.0 

2.8 

E-4 

10 Feb 68 

3.9 

2.45 

1.09 

0500 

0.9 

2.3 

E-10 

10 Feb 68 

1.4 

4.3 

11.6 

1300 

0.3 

3.3 


12 Feb 68 

1.9 

0.94 

0.65 

1100 

0.6 

1.85 

E-ll 

10 Jun 68 

1.2 

525. 

398. 

1230 

0.3 

501. 


10 Jul 68 

0.8 

141. 

252. 

0300 

0.8 

112. 


2 Sep 68 

0.6 

3.84 

1.26 

1200 

1.0 

3.02 


28 Sep 68 

0.8 

67. 

57. 

0800 

0.8 

75. 


28 Sep 68 

0.6 

162. 

119. 

2000 

1.0 

163. 


30 Sep 68 

1.2 

350. 

305. 

2100 

0.3 

316. 


30 Sep 68 

1.2 

285. 

213. 

1400 

0.4 

282. 


30 Sep 68 

0.8 

390. 

253. 

1600 

0.8 

385. 
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observed prior to a reversal in the interplanetary magnetic field are 
not present in the data within 0(1 hour) after the reversal. Additionally, 
there appears to be no distinction in this regard between persistent HPL 
depressions and persistent HPL enhancements. 

This strong correlation between the interplanetary sector and the 
persistent features, embodied in the a-pole/B-pole terminology, can now 
be used to help organize and interpret additional aspects of the proton 
observati ons . 

EDP Events 

The observations of EDP events are of particular interest because 
these events can be associated with an interplanetary region of limited 
spatial extent which is propagating radially away from the sun (and thus 
down the length of the geomagnetic tail) at a reasonably constant and 
well defined velocity: that of the solar wind. As a consequence, time 
delays between the appearance of the event in different polar cap regions 
must be directly related to the time required for the region of enhanced 
flux in interplanetary space to propagate from one access window to the 
next. 

A total of eleven EDP events and possible EDP events have been ob- 
served with this experiment. These are tabulated in table VII-3, which 
includes an indication of some of the factors leading to the identification 
of each event as an EDP event. In every case the event was observed first 
at low polar latitudes, which recent work [26,95] involving the 
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TABLE VI 1-3 

0G0-4 Observations of EDP Events 


Date 

Univ. 
Time 
of LPL 
Peak 
(HHMM) 

11 Aug 67 

0400 

19 Aug 67 

0548 

10 Nov 67 

2145 

15 Nov 67 

0400 

15 Nov 67 

0830 

1 Dec 67 

2130 

2 Dec 67 

1300 

17 Dec 67 

0600 

30 Dec 67 

1135 

30 Dec 67 

1420 

27 Aug 68 

1515 

1 Oct 68 

0130 


Average delays: 


Delay (hrs. UT) 
between LPL 
Peak and: 


a-HPL 

Peak 

8-HPL 

Peak 

0.0 

2.0 

1.2 

13. 9° 

1.1 

—d 

e 

6.0 

e 

[ 

5.1 

1.2 ■ 

7.0 

1.3 C 

3.8° 

0.0 

21.2 

1.4 

5.9 

2.2 

e 

0.4 

1.3 


5.0 


1.4±0.4 6.9±3.3 


mrt SC/SI 


Neutron 
Monitor 
Indi- 
cati on? a 

within 
4 hrs. 
of LPL 
Peak? 

FD 

YES 

WD 

YES 

NO 

NO 

NO 

NO 

NO 

NO 

WD 

NO 

NO 

NO 

NO 

NO 

FD 

YES 

FD 

YES 

NO 

NO 

FD 

YES 


Unam- 
biguous 
I dent, 
as EDP 
Event?^ 

Figure 

Showing 

Profile 

NO^ 

A-l 

YES 

A- 2 

YES 

A- 4 

YES 


NO 


YES 

A- 5 

YES 

A-6 

NO 


YES 

A- 7 

YES 

A- 7 

NO 


YES 



a. FD = Forbush Decrease; WD = weak depression 

b. All events in this table exhibit most of the phenomenological 
characteristics of EDP events (e.g. only observed at low energies, 
rapid variations in counting rate, short duration), but those indi- 
cated by "NO" in this column have some facet which is not completely 
consistent with EDP observations {Cf. discussion in Appendix A con- 
cerning the 11 August 1967 event). 

a. Observation is uncertain or impossible due to a data gap. 

d. Sector reversal occurs before 8-HPL EDP peak would be expected. 

e. Observation impossible due to orbit which did not penetrate to HPL. 
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calculation of proton trajectories in a model magnetosphere has indicated 
corresponds to an access window in the vicinity of the earth (within about 
20 After a short delay, which averaged 1.4 ± 0.4 hours, the event 

W 

was observed at a-high polar latitudes. After a much longer time delay 
(6.9 ± 3.3 hours) the event was observed at e-high polar latitudes. The 
following conclusions can be made and are illustrated in figure VI 1-2 (an 
average solar wind velocity of 400 km/sec is assumed in calculating the 
access window "positions"): 

1. THE ACCESS WINDOW FOR 1.2-40 MeV PROTONS OBSERVED 
AT a-HIGH POLAR LATITUDES IS SOMEWHAT FURTHER FROM THE SUN 
THAN THE ACCESS WINDOW FOR THOSE PROTONS OBSERVED AT LOW 

POLAR LATITUDES, (a-high polar latitude access window is located 
approximately 320 ± 90 R^ behind the earth.) 

2. THE ACCESS WINDOW FOR 1.2-40 MeV PROTONS OBSERVED 

AT 6-HIGH POLAR LATITUDES IS SIGNIFICANTLY FURTHER FROM ' 

THE SUN THAN THE ACCESS WINDOW FOR THOSE PROTONS OBSERVED 
AT LOW POLAR LATITUDES. (6-high polar latitude access window 
is located approximately 1560 ± 750 Rq behind the earth.) 

Examples of these EDP observations are given in figures A?1 , A-2, A-4 
through A-7, and V-4. 

In addition to the location of the access windows, estimates of the 
extent of the windows can also be inferred from these data. Although in 
general the polar cap flux is related to the interplanetary flux by a 



Figure VI 1-2 


Schematic representation of the configuration of the magnetospheri c 
access windows for 1.2-40 MeV protons. This configuration is referred 
to as the Lag configuration in the text. This diagram is intended to 
be completely model -independent; magnetic fields are not indicated, nor 
are details concerning the method by which protons gain access. 
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convolution integral, all that is needed at this point is an estimate. 

For this purpose, the "width" of the feature, AT, as shown on a profile, 
may be related to the spatial extent of the feature in interplanetary 
space, L ip , and the spatial extent of the access window, L w , by the approx- 

i mati on 



* V e 9 = at 
v sw 


reg 


(7.1) 


where V sw is the solar wind velocity, and "reg" represents one of "LPL", 
"a-HPL", or "e-HPl". Using M hr. as a typical value for L ip (indicated 
by the data presented by McCracken, et al- [66] and Anderson [47]), the 
0G0-4 EDP observations give the following estimates for the extent of the 
access windows: 


L LPL % 50 R ® 

K.-HPL ' 110 R ® (7 ' 2) 

L g-HPL % 140 R ® 

These numbers are very rough estimates due to the sparcity of well defined 
data and the assumptions made above, but the ratio between any two of 
these values can be considered accurate to within a factor of 2. 


Questions relating to the details of the magnetic field configuration 
in the vicinity of the access windows and the method by which this access 
is achieved will be discussed below, after considering the extent to which 
the solar flare event observations support the Lag access window config- 
uration illustrated in figure VI 1-2. 
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Fleece Events 

The interpretation of the EDP event observations was aided in no 
small degree by the availability of several events for which the data 
coverage was virtually continuous and almost uninterrupted. Because of 
the much longer duration of solar flare particle events, conclusions drawn 
from observations of flare events are, in contrast, subject to the fol- 
lowing two limitations: (1) the time-sharing nature of the 0G0-4 tele- 
metry made it virtually impossible to obtain a profile of a complete 
flare event (see Section III), and (2) the interplanetary configuration 
of proton fluxes during solar flare events is highly dependent on non- 
uniform interplanetary parameters and parameters related to the parent 
flare. The first point implies that several flare events must be ob- 
served to construct a comprehensive picture of particle access. The 
second point, on the other hand, implies that this is a limitation not 
confined to data from 0G0-4. Indeed, any conclusions drawn from obser- 
vations of a single flare event are subject to the limitation that only 
one set of solar flare and interplanetary parameters are represented. As 
the basis for the conclusions presented here, though, more than twenty- 
five solar flare events were investigated, representing a wide range of 
these parameters; these observations are summarized in table A-l. 

Although interpretations of the proton flare event observations are 
complicated by the effects of interplanetary propagation, at least four 
aspects of these observations imply an Lag configuration for the access 
regions: 
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1. Flare onsets consistently occur later at e-high polar 
latitudes than at either low polar latitudes or a-high polar lat- 
itudes. The 2 November 1967 event (figures V-7 and A-3) is a par- 
ticularly good example of this: there are two e-pole observations 
after the beginning of the flare is seen at LPL, showing a flux com- 
parable to the pre-event flux at e-HPL. The 13 July 1968 event, 
shown in figure A-12, is another example of this delayed onset 
phenomenon. 

2. Flare onsets are sometimes observed to occur later at 
a -HPL than at LPL. Once again, the 2 November 1967 event (figures 
V-7 and A-3) is a good example of this type of observation. 

3. The "shape" of the g-HPL profile relative to the LPL 
profile consistently exhibits the characteristics expected in obser- 
vations taken further from the sun: in all cases the g-HPL peak (a) 
is broader, (b) represents a smaller maximum flux, and (c) occurs at 
a later time. The 13 July 1968 event (figure A-12) is a perfect ex- 
ample of this phenomenon. Note that the e-HPL peak is clearly broader, 
lower, and later that the LPL peak. Other complete examples are 
scarce due to the data acquisition problem, but the 1 February 1968 
event (figure A-8) is notable. 

4. Occasionally the peak e-HPL flux is greater than the LPL 
flux measured at the same time. This results in the persistent feature 
being an enhancement, normally during the decay phase of the event. 

The 1 February 1968 (figure A-8) and 13 July 1968 (figure A-12) events 
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are graphic examples of this type of observation, but the 13 May 
1968 event (figure A-ll) is particularly notable: a significant 
g-HPL enhancement lasting for more than forty hours. An access 
window configuration of the Lag type discussed above would mean 
that a g-HPL enhancement should be observed whenever the inter- 
planetary flux gradient is positive and large enough: 

1 tk- > 5 x 10" 3 %/FL for a VjVj rate of ^10 cts/sec. 

J dK w 

The consistency with which all four of these phenomena are observed is 
strong evidence for the existence of the Lap access window configuration 
illustrated in figure VI 1-2. 

Thus, the 0G0-4 proton observations consistently imply that the 
1.2-40 MeV proton access windows are often related in the manner illustra- 
ted in figure VI 1-2, which we are referring to as the Lag configuration. 
This configuration and other aspects of these data are quite relevant to 
the problem of placing constraints on the possible configurations of the 
distant geomagnetic tail. 

Models of the Distant Geomagnetic Field Configuration 

As mentioned above, one of the important and currently unresolved 
questions concerning the geomagnetic field is the degree to which the solar 
and terrestrial magnetic fields merge. Although charged particle obser- 
vations cannot answer this question directly, a great deal of insight can 
be gained by analyzing the implications of such observations in the context 
of possible configurations of the geomagnetic field. Several such 
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configurations can be postulated, of course, but we will deal here with 
three fundamentally different models: a closed configuration (i.e., no 
merging), referred to as model A, an open configuration involving merging 
near the sub-solar point on the front of the magnetosphere ( model B) , and 
an open configuration involving merging at the polar neutral points ( model 
c). Each of these models will be dealt with in turn: after a brief descrip- 
tion of the major facets of a model, some of the consequences with respect 
to charged particle observations in the polar caps will be investigated. 
Following the discussion of the model, the pertinent 0G0-4 data will be 
compared to these predictions. The results of this comparison between 
the model predictions and the 0G0-4 results are listed in table VI 1-4. 

In this table, the predictions which are not attributable to one or more 
of the proponents of a given model and are thus the responsibility of 
this paper are indicated by italics; these extrapolations will be dis- 
cussed below. 

Closed Magneto spheric Configuration — MODEL A 

In the absence of merging, there will be no direct link between the 
geomagnetic field and the interplanetary magnetic field. Michel and 
Dessler have dealt with this picture extensively [3,10,96-100]. The 
most striking characteristic of this model is a very long tail, extending 

I 

as far as an astronomical unit or more behind the earth [3]. As one moves 
away from the earth along the tail, the tail is pictured as being gradually 
flattened by the anisotropic pressure of the interplanetary magnetic field 
looping over the tail. It is postulated that somewhere between 1.1 and 
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1.5 AU the tail has become so drastically flattened ('vlO 3 wide by 
'vl 0 R| thick at about 10 4 R^ behind the earth) that the field becomes 
filamentary [10]. It is questionable, however, whether such a degree of 
flattening is supportable [101]. 

The primary mode of access for charged particles is by diffusion 
through the sides of the tail [10,97-99]. Michel and Dessler also pos- 
tulate that charged particles may be able to rapidly gain access to the 
geomagnetic field in the filaments (i.e., beyond 1.1-1. 5 AU) [10]. Michel 
and Dessler further point out that such a mechanism will probably be more 
effective for one pole than the other, depending on the orientation of the 
interplanetary field near the pertinent filaments in such a way that par- 
ticles should be able to enter the filaments of the a-tail (defined in 
terms of the sector at 1.1-1. 5 AU) more readily than those of the p-tail. 

This could lead to north/south differences in the fluxes observed at high 
polar latitudes, which would be correlated with the configuration of the 
interplanetary field at 1.1 -1.5 AU [cf. A1 and A3 in table VI I -4). The 
north/south component of the field does not appear to be pertinent to these 
access mechanisms ( Cf . A2). Although these mechanisms appear to be inde- 
pendent of the presence of interplanetary anisotropies ( Cf . A8, A9 and A10), 
it would seem that interplanetary flux gradients between 1 AU (LPL access 
window) and 1.1-1. 5 AU could result in the observation of persistent features 
in the polar caps {cf. A7). 

It is instructive to cast these access mechanisms in the framework of 
access windows. Both poles are said to have a window for slow, diffusive 
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access extending from £1 All to 51.1-1.5 AU, in which the rate of access 
is postulated to be independent of the orientation of the interplanetary 
field, and hence the same in both windows. Beyond these access windows, 
in the filamentary portion of the tail ( Cf . All), is the rapid access 
window, which is postulated to be effective only for the a-pole. This 
access window configuration implies that the 3-HPL flux can never exceed 
the a-HPL flux, since any flux gaining access to the 3-tail would gain 
access to the <x-tail at the same rate through the diffusive windows [cf. 

A6) . The predictions of this model with respect to the observation of 
EDP events also follow from this picture of the access windows. The EDP 
peak would, of course, be observed at LPL coincident with the arrival at 
1 AU of the interplanetary region of enhanced flux. After a delay of 
10-50 hours (approximate solar wind propagation time for 0.1 -0.5 AU), 
a peak would be observed at a-HPL. If the diffusion rate in the dif- 
fusive access windows were sufficiently rapid, an additional a-HPL EDP 
peak might be observed, coincident with a 3-HPL peak and before the a-HPL 
peak due to access through the rapid access window. Whether or not the 
diffusion rate is rapid, however, it is clear that in either case the EDP 
observations would be phenomenologically quite different than those ob- 
served with 0G0-4 ( Cf . A4 and A5). 

As indicated in table VI 1-4 these predictions are inconsistent with 
the 0G0-4 data. With the severe constraints which these data place on this 
model, it is questionable whether it can remain supportable in its present 
form. This does not imply, however, than no other closed magnetospheric 
configuration could be developed which would meet these constraints. 
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Magnetic Merging at the Sub-solar Point — MODEL B 

In contrast to the situation in model A, significant merging between 
the geomagnetic field and the interplanetary magnetic field will result 
in a configuration in which the high polar latitude regions are directly 
connected to interplanetary space. In 1961, Dungey [2] pointed out that 
in the presence of a southward component in the interplanetary magnetic 
field there might be two magnetic neutral points in the geomagnetic- 
interplanetary magnetic field system: one on the solar side of the mag- 

netosphere, perhaps near the sub-solar point on the magnetopause, and the 
other behind the earth. As he and others have pointed out [102-107], if 
conditions at the neutral points were such as to allow magnetic field 
merging, then a configuration could be postulated wherein geomagnetic field 
lines merge with interplanetary field lines on the solar side of the mag- 
netosphere, are pulled back across the polar caps as the solar field con- 
tinues to propagate with the solar wind, and eventually reconnect at the 
anti -solar neutral point, after which the lines migrate back around the 
polar caps to the front of the magnetosphere again. Since the length of 
the tail in this configuration would be the product of the "age" of the 
tail (average time from merging at the neutral point on the solar side of 
the magnetopause to reconnection at the anti -solar neutral point) and the 
average propagation velocity, near the magnetopause, of the plasma in which 
the interplanetary field line is imbedded, the north and south geomagnetic 
tails would necessarily be the same length. Although both are the same 
length, as Dungey [106], Reid and Sauer [13], Van Allen, et at. [25], and 
others have pointed out, one tail may be connected to portions of the 
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interplanetary field which are solar-directed, while the other is con- 
nected to an anti -solar-directed field. This effect depends on the orien- 
tation of the solar field in the ecliptic plane (e.g. the sector of the 
solar field) near the earth (cf. B1 and B3 in table VII-4), and the entire 
process is dependent on the presence of a southward component in the solar 
field (Cf. B2) . 

The access mechanism for low rigidity particles in this model is 
rather simple: the direct connection between the polar cap magnetic field 
and the interplanetary magnetic field allows interplanetary fluxes to gain 
access to the polar cap regions directly, although the trajectories in- 
volved may be non-adiabatic. Reid and Sauer [13], Van Allen, et at - [25], 
and others [24,26,108] have pointed out that interplanetary anisotropies 

may contribute to features observed in the polar cap regions (cf. B8), 
especially if adiabatic motion is invoked to imply that the only particles 
seen at a-HPL will be those having interplanetary pitch angle within a few 
degrees of zero, while the only particles seen at b-HPL will be those having 
interplanetary pitch angles within a few degrees of 180°. The validity of 
this assumption is discussed below and in Appendix B. 

Several consequences for polar cap charged particle observations can 
be extrapolated from this model; most of these are consequences of the im- 
plied access window configuration: both HPL access windows located the same 
distance (the most recent estimates being a few hundred earth radii [25]) 
behind the LPL access window, which is, again, within ^20 of the earth 
( Cf . Bll). Although interplanetary flux gradients could not, of course, 
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cause north/south asymmetries, they might cause persistent features if 
large enough J- ^ > 3 * lO" 2 to cause discernible differences 

between the LPL flux and the HPL flux ( Cf . B7). Clearly, then, with that 
one possible exception, persistent features and north/south asymmetries 
would not be expected except in the presence of interplanetary anisotro- 
pies (Cf. B9 and BIO). Also, if the interplanetary anisotropy were solar- 
directed, one would expect to observe a larger flux at 6-HPL than at 
a-HPL (Cf. B6) . This access window configuration also implies that dur- 
ing EDP events both HPL peaks would be observed at the same time: on the 
order of an hour after the LPL peak ( Cf . B4 and B5). 

The disappearance of persistent features may or may not be correlated 
to sector reversals, depending on the effect a given reversal has on the 
magnitude of the interplanetary anisotropy. If the anisotropy disappears 
coincident with the sector reversal, then so also will the persistent fea- 
ture. If the anisotropy is unaffected by the sector reversal, then the 
feature should be seen to change from one pole to the other, staying in the 
6-pole (or cx-pole) ( Cf . B3). 

The extent to which the 0G0-4 data agree with these predictions is 
illustrated in table VII-4. The net result of the observations is the im- 
plication that although interplanetary anisotropies may contribute to the 
features observed at the polar caps, other causes are clearly important as 
well. In particular, the EDP event observations are somewhat difficult to 
interpret in the context of an open field configuration in which the north 
and south HPL access windows are the same distance from the earth. 
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Magnetic Merging at the Polar Neutral Points — MODEL C 

Frank has recently proposed a new variation of the open geomagnetic 
field configuration in connection with his studies of the relationship 
among the polar cusp region, the neutral sheet region, and auroras [4,109, 
110]. Although this model was originally intended to explain other phen- 
omena, its relevance to the problem of particle access into the magneto- 
sphere is unique: it is an open configuration in which widely separated 
access windows might be a natural consequence. It is important, therefore, 
to investigate what constraints must be placed on the model to yield the 
sector-dependent Lag access window configuration prior to discussing the 
consequences of the model with respect to charged particle access. First, 
though, a brief description of the model is necessary. 

All merging which takes place between the geomagnetic field and the 
interplanetary field is assumed to take place only at the polar neutral 
points (see figures VI-3 and VII-3). After a geomagnetic line of force 
merges with an interplanetary line of force, the foot of the line of force 
is pulled back along the boundary of the high polar latitude region. The 
intersection of the line of force and the magnetopause travels along an ex- 
tension of the polar "cusp" region along the flank of the magnetosphere 
until the neutral sheet region is reached. At this point one of two things 
occurs, depending on the pre-merging configuration of the geomagnetic field 
line: 

A. The geomagnetic field line was originally a closed field 
line. The other end of the original line, which connects to the other 
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Figure VII-3 

Schematic representation of the "history" of an open geomagnetic field 
line which merges with an interplanetary magnetic field line at the 
northern polar neutral point. The geomagnetic line of force merges at 
the neutral point (a), is pulled back along the cusp region (b) and 
into the neutral sheet region (c). The line eventually begins to migrate 
across the polar cap (d), until it once again passes near the polar 
neutral point (e). 
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pole, has been following a similar path along the "cusp" region in 
the other hemisphere. In the neutral sheet the lines again merge, re- 
connecting the two ends of the original interplanetary field line. The 
reconnected geomagnetic field line is eventually convected back to the 
front of the magnetosphere to undergo the process again. 

B. The geomagnetic field line was originally an open field 
line. The other end of the original geomagnetic field line is no 
longer connected to the earth, and consequently appears to play no 
further part in the process of particle access. The line of force con- 
nected to the earth eventually begins to migrate across the tail until 
it is passing near the neutral point again, at which time the pro- 
cess can start again. The migration takes place to replenish open 
field lines which are being removed from the neutral point by the 
merging process (see figure VII-3). 

Since only the open field lines are involved in the definition of access 
into the high polar latitude regions, we will confine our attention to the 
second process. 

Since the interplanetary lines of force continue to be convected away 
from the earth by the solar wind, the length of the geomagnetic tail would 
be proportional to its "age": i.e., the time required for a line to complete 
an entire cycle, from merging to merging again. This time, in turn, is in 
some sense inversely proportional to the rate at which open field lines 
merge at the neutral point. It is immediately obvious, of course, that if 
the open field line merging rate were different at the two neutral points. 
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then the north and south geomagnetic tails would be of different lengths. 
There is, therefore, a mechanism whereby an Lag configuration could be 
established. It remains, however, to investigate what conditions must be 
met to achieve the configuration of figure VI 1-2, which the 0G0-4 data have 
established. There are at least three features which are pertinent in this 
regard: 

(1) the relative locations of the access windows, 

(2) the extent of the access windows, and 

(3) the time correlation between the disappearance of persistent 
features and sector reversals in the interplanetary field. 

(1). The observational evidence which was discussed in relation to 
the Lag access window configuration illustrated in figure VI 1-2 indicates 
that the position of the B-HPL access window can be between 3 and 10 times 
as far behind the earth as the a-HPL access window. For this to occur, 
the rate at which open field lines merge at the a-polar neutral point must 
be 3-10 times the rate at the B-polar neutral point. The theoretical in- 
vestigation of this point is rather complex and is considered in detail 
in Appendix C. One of the major difficulties encountered in trying to 
deal with magnetic merging in this configuration is the applicability of 
earlier studies, all of which have concentrated on a configuration like 
that shown in figure VI-4: two exactly antiparallel fields considered in 
a plane parallel to the fields and perpendicular to the interface between 
them. The work of Sweet [84], Petschek [86], Sonnerup [87], and Yeh and 
Axford [88], discussed above, can be generalized to treat two non- 
antiparallel fields by the superposition of a constant field perpendicular 
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to the plane in which they analyze the antiparallel field configuration 
and by making the suitable minor readjustments of their results (see 
(6.1) and (6.2)). The configuration of the fields at the geomagnetic 
neutral points, however, involves, as shown in figure C-2, the interac- 
tion between a region where the field is homogeneous and a region where 
the orientation of the field is a strong function of position. Assump- 
tions must consequently be made concerning the merging likelihood for 
two fields as a function of the angle between them in the presence of 
fields at all such angles. The implications of three such assumptions, 
chosen as representative of a wide range of possible assumptions, are 
investigated in Appendix C. 

The main conclusions to be drawn from the results of Appendix C are 
that (a) such differences in polar merging rates are well within the limits 
of reasonable parameters if certain assumptions are made concerning the 
angular dependence of the merging rate, and (b) if the sector is redefined 
in terms of the earth-sun line and the interplanetary field direction near 
the polar neutral point instead of the average spiral angle and the field 
direction far from the earth, the correct sector correlation follows. 

That the critical interplanetary field configuration is that of the 
portion of the field near the polar neutral points complicates matters con- 
siderably. It is not clear how this configuration depends in detail on the 
field configuration in interplanetary space and on the solar wind parameters 
[Cf. 111,112]. Nor is it clear that for a given set of interplanetary 
parameters the interplanetary field will make the same angle with the 
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earth-sun line at one neutral point as it does at the other. Until this 
question is resolved, it will not be possible to determine whether there 
is an exact, detailed correlation between the access window positions 
inferred from the data and those inferred from this model. 

(2). The size of the HPL access windows relative to their positions 
places constraints on the way in which the geomagnetic field lines mi- 
grate from the neutral point to the nuetral sheet region and back to the 
neutral point. Figure VI I -4 shows a tentative migration pattern, wherein 
the paths followed by the feet of the open lines of force are shown pro- 
jected into an invariant latitude-magnetic local time coordinant system. 
After merging at the neutral point a line is pulled down the cusp region 
(paths A), into the neutral sheet region (paths B), and finally migrates 
back to the neutral point (one of paths C). As soon as the foot of a 
line of force enters the central region of the polar cap (C), the inter- 
section of that line with the magnetopause becomes part of the HPL access 
window for that pole. It remains part of the access region until the 
magnetopause intersection is no longer connected to the polar region 
(i.e., until the line merges again at the neutral point). The length of 
the access window (measured parallel to the axis of the tail) is there- 
fore proportional to the time required for the lines to migrate from the 
neutral sheet region (B) across the polar cap (C) to the neutral point. 

The data presented in Section V indicate the parameters given in 
(7.2) as "typical" for the Lag access window configuration. In contrast 
to the case for the positions of the access windows, the length of the 
windows does not seem to depend heavily on the interplanetary sector. This 
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Figure VI 1-4 

Tentative sketch of the migration pattern for polar geomagnetic 
field lines in Frank's magnetospheric model (model C). Examples 
of the "paths" followed by the feet of the polar field lines are 
projected into a A-MLT coordinate system. 
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would imply that the time spent by a line migrating across region C would 
need to be more or less independent of the total migration time for a 
complete cycle. This is a clear constraint which the experimental data 
place on this configuration. 

(3). The data shown in table VII-2 indicate that persistent B-HPL 
features disappear within about one hour after a reversal in the inter- 
planetary sector. This means that within one hour after a sector reversal 
from, for instance, negative to positive the northern high polar latitude 
access region has changed from a e-HPL configuration to an a-HPL config- 
uration. This establishes a lower limit for the merging rate at the 
a-polar neutral point: it must be at least as large as that necessary for 
almost every high polar latitude field line to have merged with the inter- 
planetary field at least once in less than one hour; i.e., an open field 
merging rate not less than ^1.3 * 10 5 uG-R|/hr. 

This rate is not inconsistent with the following estimate for the 

i 

maximum rate at which the interplanetary magnetic field impinges on the 
polar neutral points. If all of the ^50 uG interplanetary magnetic field 
which impinges on the front of the magnetosphere is thought of as being 
swept up past the neutral point by the ^225 R # /hr (400 km/sec) solar wind, 
then the "rate" at which the interplanetary field would impinge on each 
neutral point would be M.4 * 10 5 yG-R|/hr. 

These rates are also consistent with at least two of the theoretical 
maximum merging rates discussed in Section VI: those predicted by Sonnerup 
[87] and by Yeh and Axford [88]. The consistency with the latter predic- 



tion is self-evident, since they predict that there will be no upper limit 
on the merging rate. In order to check the consistency of these merging 
rates with the prediction of Sonnerup, given by (6.2), we must estimate 
the solar wind conditions near the polar neutral points. According to 
Spreiter, et aZ. [Ill], 

V _/V * 0.7 

np sw 

and 

p__/p_ tl ^ 1 .8 

np sw 

Also, according to Alksne [113], 


This means that 
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and the Alfvenic Mach number of the solar wind plasma near the polar 
neutral point is therefore 


<Vnp ■ <VswV 


V /V 
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— nfr — 
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^ 1.9 


Thus the idea of the interplanetary magnetic field merging as fast as it 
impinges on the neutral point is not inconsistent with (6.2), which specifies 
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Assuming these constraints to be satisfied, the degree to which the 
0G0-4 data agree with the resultant predictions of this model is indicated 
in table VII-4, although the a priori assumption that the constraints dis- 
cussed above were satisfied was begging the question somewhat. It is inter 
esting that this model is the only type discussed here which even quali- 
tatively provides a mechanism for generating an Lag access window config- 
uration. Although the data presented here cannot provide a conclusive 
test of the validity of the model in its present form, they do provide 
three major constraints within which more refined versions of the model 
will have to be developed. 

Interplanetary Anisotropies 

Although not directly supportable from the 0G0-4 data due to a lack 
of interplanetary measurements, it is certainly not unreasonable to suppose 
that interplanetary anisotropies may have some affect on polar cap obser- 
vations and may contribute to persistent features ( Cf . #8 in table VII-4). 
The entire role of interplanetary anisotropies in the formation of polar 
cap features in an open geomagnetic configuration is not at all as straight 
forward as it may appear at first, however. The assumption of adiabatic 
motion in the propagation of particles from interplanetary space to the 
polar regions means that the magnetic moment is conserved and that there- 
fore the only particles observable at the altitude of 0G0-4 are those with 
interplanetary pitch angles sl° or 5:179°. Since the geomagnetic tail 
field and the interplanetary field are almost never parallel, the propa- 
gation of particles from interplanetary space to the geomagnetic tail 
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will involve the passage through a region where the magnetic field is 
changing direction. Since the propagation of particles in the tail be- 
tween the transition region (access window) and the polar caps can be 
considered adiabatic, only those particles with pitch angles near zero 
(or 180°) in the tail will be observed at the altitude of 0G0-4. If the 
"curvature" of the magnetic field in the transition region is sufficiently 
large compared to the gyroradius of the particle, then adiabatic motion 
might indeed be expected. 

Appendix B treats this problem in some detail, dealing with the con- 
figuration of two magnetic fields equal in magnitude separated by a re- 
gion in which the field changes direction with a constant radius of 
curvature {Cf. figure B-l). Figure VII-5 shows typical results for the 
magnetic fields at an angle of 60°. In this figure, <j> ip is the inter- 
planetary pitch angle giving a pitch angle of 0° after traversing the 
transition region; <j> ip is shown as a function of k/p, where p is the 
particle's gyroradius and < is the constant radius of curvature of the 
field in the transition region. Other values of y give comparable results 
(see Appendix B). These results show that, for y=60°, unless k/p > 50, 
the interplanetary pitch angle ultimately being observed at the polar 
caps varies drastically as a function of particle energy and as a func- 
tion of the instantaneous interplanetary field configuration. One would 
not expect, for instance, that 3 MeV protons, which have a gyroradius 
of about 9.8 R^ in a 60 yG interplanetary field, with pitch angles <1° in 
interplanetary space would be observed at the polar caps, except transi- 
torily, unless k were greater than about 490 R ffi . Since the cross sectional 
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Figure VII-5 

Interplanetary pitch angle (.$^ ) giving a pitch angle of 0° after 
passing through a region in which the magnetic field changes direction 
through an angle y=60°. The direction change is assumed to take place 
at a constant radius of curvature, k . The gyro-radius of the particle 
in the magnetic field, whose magnitude does not change, is represented 
by p. See Appendix B for a more detailed explanation of these results. 
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radius of the entire tail at about 200 is about 25 Rg, it is hard to 
visualize a transition region which is at least 490 R^ wide. These pitch 
angle mapping results are, of course, relevant only in the presence of an 
interplanetary anisotropy; if the interplanetary flux is isotropic, it does 
not matter which interplanetary pitch angles are observed over the polar 
caps . 


Although direct observational evidence for interplanetary anisotropies 
is not available from the 0G0-4 data, there are three pieces of pertinent 
indirect evidence indicating the observation of persistent north/south 
differences in the absence of interplanetary anisotropies ( Cf . #10 in 
table VII-4): 

(a) : A number of the flare events observed with 0G0-4 either were associ- 

ated with east limb optical flares or were possibly associated with 
optical flares on the other side of the solar disc. Although inter- 
planetary anisotropies have been observed at low proton energies during 
same east limb flares, it is suggestive that persistent features (and 
hence north/south flux differences) were observed during all flare 
event observations made with 0G0-4. 

(b) : Persistent features are consistently observed during the decay phase 

of flare events, whereas anti -solar directed field aligned inter- 
planetary anisotropies normally disappear shortly after the peak of 
a flare event. Interplanetary anisotropies observed during the decay 
phase of flare events are normally anti -solar directed and solar wind 
aligned. Examples of persistent e-HPL depressions are shown in fig- 
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ures A-3 and A-9. The g-HPL flux was as much as a factor of 10 
below the LPL and a-HPL fluxes during the 8 February 1968 event 
(figure A-9). 

(c): Several features are observed to persist for very long periods, 

exceeding 40 hours in at least two of the cases shown in Appendix A. 
This is a much longer duration than that expected for large interplan- 
etary anisotropies (5-15 hours, depending on the time between onset 
and the peak). Figure A-3 shows a persistent feature which lasted for 
more than 24 hours. The best example, though is probably that in fig- 
ure A-9, showing a very large feature ((V^)^: (Vi\Q b _ hpl as 
large as 25:1) persisting for at least 31 hours during the decay 
phase of the flare event. 
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VIII. CONCLUSIONS 


Observations of polar cap fluxes of low energy solar protons 
(1.2-40 MeV) and electrons (0. 4-1.0 MeV) made with an experiment on 
board the 0G0-4 satellite lead to the following conclusions. 

1. Extensive observations of electron polar cap fluxes during 
geomagnetically quiet periods have resulted in a comprehensive mapping 
of the boundary between open and closed geomagnetic field lines onto 
the polar cap. This mapping, shown in figure V- 3, is an order of mag- 
nitude more comprehensive than previously available. 

2. Observations during 54 solar proton events have established 
that the presence of persistent features in the polar cap proton flux 
is the norm rather than the exception (see table A-l and figure VII-1). 

3. A comparison between the electron polar cap boundary (see con- 
clusion #1) and persistent features in the proton polar cap flux indi- 
cate that, in general, these features represent different flux levels in 
the open field line region (high polar latitudes — HPL) than in the 
closed field line region (low polar latitudes -- LPL). 

4. There is a strong correlation between the observation of persis- 
tent polar cap features and the sector structure of the interplanetary 
magnetic field: in general, persistent features are observed in the south 



150 


pole during a positive sector and in the north pole during a negative 
sector. On the basis of this correlation, the term "3-pole" is defined 
in such a manner as to refer to the south pole during a positive sector 
and in the north pole during a negative sector; in both cases the op- 
posite pole is referred to as the a-pole. The 0G0-4 observations reported 
here are correlated in this fashion 91% of the time. 

5. The observation that persistent features disappear within o(one 
hour) of a reversal in the interplanetary sector is also consistent with 
this sector correlation. 

6. There appears to be little or no correlation between the presence 
of persistent features in the proton polar cap flux and the orientation of 
the north/south component of the interplanetary magnetic field: during the 
observations of persistent features this component was southward only 54% 
of the time and changed orientation several times with no noticeable 
effect. 

7. Observations of solar proton events, especially those associated 
with co-rotating regions of enhanced flux in interplanetary space, have 
been used to establish the La3 access window configuration illustrated in 
figure VI 1-2, in which the LPL access window is near the earth, the a-HPL 
access window is a few hundred earth radii behind the earth (in the anti- 
solar direction), and the s-HPL access window is located <9(2000 Rj behind 
the earth. 

8. The widths of the HPL peaks observed on the EDP event profiles 
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imply that the extent of the a-HPL access window is about the same as 
that of the e-HPL access window. A comparison with the LPL peaks further 
implies that neither of the HPL access windows is much larger than the 
LPL access window (see figure VI 1-2). 

9. Although direct interplanetary measurements are not available 
from the 0G0-4 data, it is concluded that the following contribute to 
persistent polar cap features (north/south asymmetries): 

A. limited interplanetary regions of enhanced flux propa- 
gating past the earth, 

B. radial gradients in the interplanetary flux, and 

C. anisotropies in the interplanetary flux. 

TO. The 0G0-4 observations are not consistent with the particle 
access mechanisms proposed by Michel and Dessler [3]. 

11. The access window configuration deduced from EDP event obser- 
vations (conclusion #7) is not consistent with an access window config- 
uration in which both HPL access windows are the same distance from the 
earth, as would be expected in an open magnetospheric configuration in- 
volving magnetic merging near the sub-solar point ( Cf . Dungey [2]). 

12. An Lag access window configuration (conclusion #7) may be a 
consequence of an open magnetospheric model involving magnetic merging 
at the polar neutral points ( Cf . Frank [4]). To achieve this configur- 
ation, however, it is necessary that the model satisfies at least the 
following constraints: 
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A. The merging rate for open geomagnetic field lines at one 
polar neutral point (a) must be ^5 times that at the other polar 
neutral point (e). This difference between the open field line 
merging rates must be related to the orientation of the interplane- 
tary magnetic field, as measured far from the neutral point, in. such 
a way as to be consistent with the observed correlation with the 
interplanetary sector structure (conclusion #4) and lack of correla- 
tion with the north/south component of the interplanetary field 
(conclusion #6). 

B. The time required for open field lines to migrate from 
the neutral sheet region across the polar cap to the vicinity of the 
neutral point must be about the same in both tails, so as to account 
for both HPL access windows being about the same size (conclusion #8). 

C. The rate at which open field lines merge with the inter- 
planetary field at the a-polar neutral point must be at least as 
large as that required for almost all of the open field lines to have 
merged in o(one hour), since 6-HPL features disappear within O(one 
hour) after a sector reversal (conclusion #5). 

In addition to these observational conclusions, the following con- 
clusions are implied by the results of the theoretical investigations pre- 
sented in Appendices B and C. 

Appendix B: The assumption of adiabatic motion for the access of 
1.2-40 MeV protons from interplanetary space to the interior of the 



153 


magnetosphere is not a good approximation. If the extent to which an 
interplanetary pitch angle distribution maps onto the polar cap region 
were to be determined from appropriate data, figure B-4 would indicate 
the constraints thus imposed on the configuration of the magnetic field 
in the vicinity of the access windows. 

Appendix C: The satisfaction of constraint A mentioned above in 

conclusion #12 is possible if the proper assumptions are made concerning 
angular dependence of the probability that two field lines will merge. 
This can be seen by comparing figures C-8, C-9 and C-10. 

Although the results outlined above are most consistent with a 
field configuration involving magnetic merging at the polar neutral 
points, the question of whether there is a significant degree of merg- 
ing between the geomagnetic and interplanetary magnetic fields cannot 
be unambiguously determined by charged particle measurements. These 
observations do, however, provide a suitable platform upon which to 
base theoretical investigations of this question. 



